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ABSTRACT 

We present the first three-frequency South Pole Telescope (SPT) cosmic microwave background 
(CMB) power spectra. The band powers presented here cover angular scales 2000 < I < 9400 in 
frequency bands centered at 95, 150, and 220 GHz. At these frequencies and angular scales, a com- 
bination of the primary CMB anisotropy, thermal and kinetic Sunyaev-Zel'dovich (SZ) effects, radio 
galaxies, and cosmic infrared background (CIB) contributes to the signal. We combine Planck/'SFl 
and SPT data at 220 GHz to constrain the amplitude and shape of the CIB power spectrum and find 
strong evidence for non-linear clustering. We explore the SZ results using a variety of cosmological 
models for the CMB and CIB anisotropies and find them to be robust with one exception: allowing for 
spatial correlations between the thermal SZ effect and CIB significantly degrades the SZ constraints. 
Neglecting this potential correlation, we find the thermal SZ power at 150 GHz and t = 3000 to be 
3.65 ± 0.69 /zK^, and set an upper limit on the kinetic SZ power to be less than 2.8 /xK^ at 95% 
confidence. When a correlation between the thermal SZ and CIB is allowed, we constrain a linear 
combination of thermal and kinetic SZ power: 

^3000 + O-S^'sTOO = 4.60 ± 0.63 AiK2, consistent with 
earlier measurements. We use the measured thermal SZ power and an analytic, thermal SZ model 
calibrated with simulations to determine (j% ~ 0.807 ± 0.016. Modeling uncertainties involving the 
astrophysics of the intracluster medium rather than the statistical uncertainty in the measured band 
powers are the dominant source of uncertainty on ag. We also place an upper limit on the kinetic 
SZ power produced by patchy reionization; a companion paper uses these limits to constrain the 
reionization history of the Universe. 

Subject headings: cosmology - cosmology: cosmic microwave background - cosmology: diffuse radiation- 
cosmology: observations - large-scale structure of universe 
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1. INTRODUCTION 

The current generation of ground-based cosmic mi- 
crowave background (CMB) experiments enables re- 
searchers to probe, for the first time, the power spectra of 
the thermal and kinetic Sunyaev-Zel'dovich (SZ) effect s 



(iSunyaev fc Zerdovich|1972| ISunyaev fc Z eldovich|1980j) 
I'he South Pole I'eiescope (SPT, Carlstrom et all 
and the A tacama Cosmology Telescope (ACi', " 



et al.|2007 ) have sufficient sensitivity, angular resolution, 
and frequency coverage to measure the power spectrum 
of the thermal SZ (tSZ) effect and to place increasingly 
strong upper limits on contributions from the kinetic SZ 
(kSZ) effect ( [Lueker et al.|2010[ |Shirokoff et al.|20l"T|[Das 



et al. 2011). 'i'hese secondary femperature anisotropies 
in the CiVit} arise when CMB photons are scattered by 
the inhomogeneous distribution of free electrons well af- 
ter recombination. As such, they depend sensitively on 
both the growth of structure and the details of reioniza- 
tion. 

The tSZ effect occurs when CMB photons inverse 
Compton scatter off hot electrons, predominantly those 
gravitationally bound in the potential wells of galaxy 
clusters. The scattering tends to increase the photon 
energy, distorting the CMB black-body spectrum with 
a loss of photons below approximately 217 GHz and an 
excess at higher frequencies. The kSZ effect occurs when 
CMB photons are Doppler shifted by the bulk velocity 
of electrons in intervening ionized gas. The kSZ effect 
changes the apparent CMB temperature while maintain- 
ing the black-body spectrum. We can discriminate be- 
tween the two effects based on their frequency depen- 
dence. 

The anisotropy power due to the tSZ effect is highly 
sensitive to the normalization of the matter power spec- 
trum, commonly parametrized by the RMS of the z = 



mass distrib ution on 8 Mpc//i scales, trg (Komatsu & 
Seljak|[2002| . As a measure of structure growth, the tSZ 



power spectrum can provide independent constraints on 
cosmological parameters and improve the precision with 
which they are determined. For instance, a measurement 
of the tSZ power spectrum can be used to determine the 
sum of the neutrino masses by breaking the degeneracy 
between erg and neutrino mass that exists with the CMB 
data alone. The tSZ power spectrum can also be used 
to constrain non-standard cosmological models, for ex- 
ample, placing hmits on the rang e of allowed early dark 
energy models (Alam et al. 2011). 



The tSZ power spectrum is challenging to model accu- 
rately because it includes significant contributions from 
galaxy clusters spanning a wide mass and redshift range. 
Modeling uncertainties arise from both non-gravitational 
heating effects in low-mass clusters and the limited body 
of observational data on high-redshift clusters. Recent 
models that vary in their treatment of cluster gas physics 
differ in amplitude by up to 50% for a given cosmolog y 



( |Shaw et al.|2010||Trac et al |2011[[Battaglia et al.|2010| 

'I'he post-reionization kSZ power spectrum is expected 
to be easier to model than the tSZ spectrum. As it is not 
weighted by the gas temperature, the kSZ power spec- 
trum depends less on non-linear effects that complicate 
models of the tSZ signal. The total kSZ power spectrum 
does, however, depend strongly on the details of reioniza- 
tion, which are not yet well understood. In the standard 



picture of reionization, ionized bubbles form around the 
first stars or quasars. These bubbles eventually overlap 
and merge, leading to a fully ionized universe. The ion- 
ized bubbles impart a Doppler shift on scattered CMB 
photons proportional to their linc-of-sight velocity. This 
so-called "patchy" reionization kSZ signal depends on 
the duration of reionization as well as the distribution of 



bubble sizes ( [Gruzinov fc Hu|1998[|Knox et al.|1998[ ). De- 
tection of the kSZ power spectrum or even an upper limit 
on its amplitude would le ad to interesting c onstraints on 
the epoch of reionization ( Zahn et al.||2005 1. 
The first detection of the combined tSZ and kSZ spec- 



tra was presented by Lueker et al. ( [2010 hereafter LIO) 
using data from the SF'l'. Since then, improved mea- 
surements of the pow er spectra have been reported by 
Shirokoff et al. (2011 hereafter Sll) based on 150 and 
220 GHz data for 200 deg^ of sky. Sll report a 4.5 ct de- 
tection of SZ power at 152 GHz at a multipole of ^ = 3000 
with an amplitude of Df^l^ -h 0.5I?|^o^o = 4.5 ± 1.0 ^K^ Q 
The 95% confidence upper limits on each component 
were Df^Q < 5.3AiK2 and D^^^^ <Q.b^iK^. The mea- 
sured tSZ power is less than predicted by some models, 
suggesting that either these models over-predict the tSZ 
power or dg is signi ficantly less than t he WMAP pre- 
ferred value of ~ 0.8 ( Larson et al.|[20lT ). 

The ACT coll aboration has al s o measured the high-£ 
power spectrum ('Das et al.|2011 ). Dunkley et al. (20111 
use the Das ct al. (2011) bandpowers to measure the sum 
of the tSZ and kSZ^wer at ^ = 3000 and 148 GHz to 
be 6.8 ± 2.9 //K2. The SPT and ACT bandpowers and 
resulting constraints on SZ power are consistent within 
the reported uncertainties. 

The sky power at millimeter wavelengths and arc- 
minute scales is dominated by the signal from bright 
synchrotron sources, with additional contributions from 
nearby infrared galaxies and a population of high- 
redshift, gravitati onally lense d, dusty, star- forming 

galaxies (DSFGs) ( [Vieira et aL[[2010[ ). After these iden- 
tifiable sources are masked, the power spectrum is domi- 
nated by the cosmic infrared background (CIB), the bulk 
of which is contributed by unlensed DSFGs below the 
detection threshold of instruments like SPT and ACT 
( Lagache et al.|2005 |. These DSFGs are important both 
as a foreground for CMB experiments and as a tool to 
understand the history of star and galaxy formation (e.g 



Bond et al.|[T986l[T99Ti|Knox et al.|[200T 
20051). The poweTIEran DSFGs has bot 



Lagache et al 
1 Poisson and 
clustered components, with distinct angular scale depen- 
dencies. Measurements of CIB power at multiple mil- 
limeter wavelengths can be combined to constrain the 
spectral index of the CIB anisotropy. Note that this de- 
pends on both the spectral energy distribution of indi- 
vidual DSFGs and their redshift distribution. Sll con- 
strained the spectral index of the DSFGs between 150 
and 220 GHz to be a = 3.58±0.09. [Dunkley et aI[pOTT | 
also detected significant power attributed to clustered 
DSFGs in the ACT data and found a preferred DSFG 
spectral index of 3.69 ±0.14, consistent with the Sll re- 



Throughout this work, the unit K refers to equivalent fluc- 
tuations in the CMB temperature, i.e., the temperature fluctua- 
tion of a 2.73 K blackbody that would be required to produce the 
same power fluctuation. "The flux conversion factor is given by the 
derivative of the blackbody spectrum, evaluated at 2.73 K. 
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suits. The Planck Collaboration et al. (2011) combined 
the Planck/ tib'l CIB bandpowers with the Sll bandpow- 
ers to show that the combined dataset is inconsistent 
with only the Poisson and linear clustering terms - non- 
linear clustering is required. 

This is the fourth SPT power spectrum analysis fo- 
cused on CIB and secondary CMB anisotro pics. Here, 
we irn prove upon the previous work by LIO, |Hall et aT 
(2010) and Sll in two ways. First, we are analyzing 
three- frequency data (with six cross-spectra) where ear- 
lier analyses used two frequency data with three cross- 
spectra. The additional frequency information allows the 
separation of CIB and SZ components. Second, we are 
analyzing four times more sky area (800 deg^) with the 
concomitant redu ction in bandpow er uncertainties. A 



companion paper, [Zahn et al.| 
derived kSZ constraints in ligi 
tion. 



Zll), interprets the 
ot the epoch of reioniza- 



The paper is organized as follows. In f|2] we present 
the analysis and the data to which it is applied. We 
examine the results of tests for systematic errors in the 
data in fjSj The resulting power spectrum is presented 
in [|4] The cosmological modeling is discussed in Sj5]and 
basic results in Spj We explore constraints on the tSZ 
effect in |7) on the kSZ effect in S|8) and on the CIB in 
^ We conclude in j flOl 



2. DATA AND ANALYSIS 

We present bandpowers from 800 deg^ of sky observed 
with the SPT at 95, 150, and 220 GHz. Bandpowers 
are estimated using a pseudo-C^ cross-spectrum method. 
We calibrate the data by comparin g the power spectrurn . 
to seven-year WMAP (WMAP7, jLarson et al]|2011 ) 
bandpowers. 

2.1. Data 

The 800 deg^ of sky analyzed in this work were ob- 
served by the SPT during the 2008 and 2009 Austral 
winters. The SPT is an off-axis Gregorian telescope with 
a 10 meter diameter primary mirror located at the South 
Pole. The receiver is equipped with 960 horn-coupled spi- 
derweb bolometers with superconducting transition edge 
sensors. The telescope and receiver are discus sed in more 
detail in[R^h l et a l. | p004l ), |Padin et"al] ( |2008[ ), and |Carl- 
strom et ar.l ( |201ip . 

Each held was observed to an approximate depth of 
18 /iK-arcmin at 150 GHz. The noise levels at 95 and 
220 GHz vary significantly between the two years as the 
focal plane was refurbished to maximize sensitivity to 
the tSZ effect before the 2009 observing season. This re- 
furbishment added science-quality 95 GHz detectors and 
removed half the 220 GHz detectors. As a result, the 2009 
data are comparatively shallower at 220 GHz, but much 
deeper at 95 GHz. The 2009 sky coverage accounts for 
approximately 75% of the total sky used in this analysis. 
When calculating the power spectra, the 2009 data has 
a statistical weight of 100% at 95 GHz, 73% at 150 GHz 
and 43% at 220 GHz. 



SPT beam at 95, 150, and 220 GHz are well-represented 
by 1.7', 1.2', and 1.0' FWHM Gaussians. As in Sll, we 
calculate the principal components of the full covariance 
matrix of this measurement and find that three eigenvec- 
tors are adequate to represent >99% of the covariance. 
We treat the eigenvectors from different years as inde- 
pendent. We have tested the opposite extreme, treating 
them as completely correlated, and have found minimal 
impact on the science results. The beam uncertainty 
is included in th e bandpower covariance matrix as de- 
scribed in |2.3.4| 

The observation-to-observation relative calibrations of 
the time-ordered data (TOD) are determined from re- 
peated ineasurementso^agalact^^ HII region, RCW38. 
As in [Keisler et al. ( 2011[ Kll), the absolute calibra- 
tion at each trequency is determined by comparing the 
SPT and WMAP7 bandpowers over the multipole range 
£ e [650, 1000]. We use the same ^-bins for each experi- 
ment: seven bins with A£ = 50. For the absolute calibra- 
tion, we do not weight indi vidual modes within each bin 
(as we will do in section 2.3 ). A uniform weighting is used 
because these large angular scales are dominated by sam- 
ple variance rather than instrumental and atmospheric 
noise. This calibration method is model- independent, 
requiring only that the CMB power in the SPT fields is 
statistically representative of the all-sky power. We es- 
timate the uncertainty in the SPT power calibration at 
95, 150, and 220 GHz to be 3.5%, 3.2%, and 4.8% respec- 
tively. These uncertainties are highly correlated because 
the main sources of error, WMAP 7 bandpower errors 
and SPT sample variance, are nearly identical between 
bands. 

2.3. Bandpower estimation 

We use a pseudo-Cf method to estimate the bandpow- 
ers (Hivon et al. 2002). In pseudo-C^ methods, bandpow- 
ers are estimated directly from the spherical harmonic 
transform of the map after correcting for effects such as 
TOD filtering, beams, and finite sky coverage. We pro- 
cess the data using a cross spectrum b ased analysis (Po- 
lenta et al.|2005 [Tristram et al.|[2005 ) to eliminate noise 
bias. Our scan strategy produces of arder 100 complete 
but noisy maps of each field, and so is ideally suited to 
a cross-spectrum analysis. We use a fiat-sky approxima- 
tion. Before Fourier transforming, each map is apodized 
by a window designed to mask point sources detected at 
greater than > 5cr (6.4 mJy) at 150 GHz, and to avoid 
sharp edges at the map borders. Beam and filtering ef- 
fects are corrected for according to the formalism in the 
MASTER algorithm ( [Hivon et aL||2002[ ). We report the 
bandpowers in terms ot 'L>i, where 



2n 



(1) 



Details of the bandpower estimator have been pre- 
sented in previous SPT papers (LIO, Sll, Kll). In the 
following sections, we will only note differences between 
the current analysis and what was done by Sll. 



2.2. Beams and Calibration 

The SPT beams are measured for each year using a 
combination of bright point sources in each field, Venus, 
and Jupiter as described in Sll. The main lobes of the 



2.3.1. Simulations 

Simulations are used to estimate the transfer function 
and sample variance. For each field, we construct simu- 
lated observations of 100 sky realizations that have been 
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smoothed by the appropriate beara. The sky reahzations 
are sampled using the SPT pointing information, filtered 
identically to the real data, and processed into maps. 

Each simulated sky is a Gaussian realization of the 
sum of the best-fit lensed WMAP7 ACDM primary CMB 
model, a kSZ model, and point source contributions. 
Note that this approach neglects non-Gaussianity in the 
kSZ and radio source contributions and ma y slightly un- 
derest imate the sample variance as a result. [Millea et al. 
(2011) argue this is negligible for SPT since the non- 
Gaussian signals are small and, on the relevant angular 
scales, the bandpower uncertainties are dominated by in- 
str umental noise. The k SZ power spectrum is taken from 
the Sehgal et al. (2010) simulations and has an ampli- 
tude of 2.05 fiK^ at £ = 3000. We include both Pois- 
son and clustered point sources. The Poisson contribu- 
tion reflects both radio source and DSFG populations. 
The amplitude of the radio source term is set by the de 
Zotti et al. (2005) model source counts to an amplituoe 



at a fixed separation from the diagonal in each block. 



D. 



3000 



1.28 ^K^ at 150 GHz with an assumed spectral 



index of a,- 



term at 150 GHz is -Dfooo 



0.6. The amplitude of the Poisson DSFG 
7.7 ^K^. Finally, the clus- 
tered DSFG component is modeled hy a D( (x £ term 
normalized to Z^foop = 5.9 /xK^ at 150 GHz. The DSFG 
terms have an assumed spectral index of 3.6. The am- 
plitude of each component was selected to be consistent 
with the Sll bandpowers. 

2.3.2. Covariance estimation and conditioning 

The bandpower covariance matrix includes two terms: 
sample variance and instrumental noise variance. The 
sample variance is estimated with the signal-only Monte 
Carlo bandpowers. The instrumental noise variance is 
calculated from the distribution of the cross-spectrum 
bandpowers jj'j^'^'^i'^^ between observations A and B, 
and frequencies Vi and Vj, as described in LIO. We expect 
some statistical uncertainty of the form 



iCtb' - {Cbb') f) = ^ 



CbbCb'b' 



nobs 



(2) 



in the estimated bandpower covariance matrix. Here, 
Hobs is the number of observations that go into the co- 
variance estimate. This uncertainty on the covariance 
is expected to be significantly higher than the true co- 
variance for almost all off-diagonal terms due to the de- 
pendence of the uncertainty on the (large) diagonal co- 
variances. We reduce the impact of this uncertainty by 
"conditioning" the covariance matrix. 

How we condition the covariance matrix is determined 
by the form we expect it to assume. The covariance ma- 
trix can be viewed as a set of thirty-six square blocks, 
with the six on-diagonal blocks corresponding to the co- 
variances of 95 X 95, 95 x 150, ... 220 x 220 GHz spec- 
tra. The bandpowers reported in Table [T] are obtained 
by first computing power spectra and covariance matri- 
ces for bins of width A£ — 200 with a total of 37 initial 
^-bins, thus each of these blocks is an 37x 37 matrix. The 
shape of the correlation matrix in each of these blocks is 
expected to be the same, as it is set by the apodization 
window. As a first step to conditioning the covariance 
matrix, we calculate the correlation matrices for the six 
on-diagonal blocks and average all off-diagonal elements 



k2—k — k' 



^kk' — 



^kl-k 



=k~k' 



1 



(3) 



We set the correlation to zero at distances £ > 400 where 
we find no evidence for correlations. We have tested that 
varying this condition by a factor of two in either direc- 
tion does not affect the cosmological constraints. This 
averaged correlation matrix is then applied to all the 
blocks. 

The covariance matrix includes an estimate of the sig- 
nal and, if both bandpowers share a common map, the 
noise variance. To illustrate the latter condition, instru- 
mental noise is included in the blocks corresponding to 
(150 X 150, 150 X 150) and (95 x 150, 150 x 220) covari- 
ances, but not the covariance between the 150 x 150 and 
220 X 220 bandpowers. Given the map filtering, we nei- 
ther expect nor observe correlated noise between different 
frequencies in the signal band. 

We impose one final condition on the estimated covari- 
ances of off-diagonal blocks for which the instrumental 
noise term is included, e.g. the (95 x 150, 150 x 220) co- 
variance. In these blocks, the noise on the covariance 
estimate for an £-hin can be large compared to the true 
covariance on small angular scales. As a result, we fit 
an approximate template based on the measured beam 
shapes from t = 2000 — 4400 and extrapolate this tem- 
plate to estimate the diagonal elements of the covariance 
matrix at ^ > 4400. We have tested zeroing these ele- 
ments of the covariance matrix instead and found it does 
not affect the cosmological constraints. 

2.3.3. Field weighting 

We construct a diagonal weight matrix for each field i, 
wli, = 1/Qf,. We smooth the weights with A£ = 1000 at 
£ < 3000 and A£ = 2600 on smaller scales to deal with 
uncertainty in the covariance estimate since we expect 
the optimal weights to vary slowly with angular scale. 

The combined bandpowers are calculated from the 
individual-field bandpowers according to: 



(4) 



where 



to 



T.<b = ^- (5) 

i 

The covariance matrix is combined likewise according 
Cbb' =^wl^Cl,,,wl,t, (6) 



2.3.4. Beam uncertainties 

The procedure for estima ting the beams and beam un- 
certainties is described in [2.2 here we discuss how the 
beam uncertainties are included in the cosmological anal- 
ysis. In Sll, three parameters per frequency band de- 
scribing the beam uncertainties were marginalized over 
in the Monte Carlo Markov chain (MCMC) with a Gaus- 
sian prior. In this work, we instead follow the treatment 
laid out in Kll and include the beam uncertainties in 
the bandpower covariance matrix. 
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Following Kll, we estimate the "beam correlation" 
matrix for each year: 



„beam 

Pbb' — 



where 



6Db 



1 



SDb 
Db 



1 



fSDb 



V Db 

SBb 
Bb 



(7) 



(8) 



and Bb is the Fourier transform of the beam map aver- 
aged across bin b. 

We combine the correlation matrix for each year ac- 
cording to the yearly weights. The correlation matrix is 
translated into a covariance matrix by multiplying by the 
measured bandpowers: 



ibcam 
'bb' 



beam 
Pbb 



DbD. 



b'- 



(9) 



This beam covariance is added to the bandpower co- 
variance matrix due to sample variance and instrumental 
noise. 

3. JACKKNIFE TESTS 

We test for systematic errors using a suite of "jack- 
knife" tests. A jackknife test consists of dividing the data 
into two halves and differencing them to remove true as- 
trophysical signals. The halves are chosen to maximize 
the signal in the differenced spectrum due to a poten- 
tial systematic bias. We then compare the differenced 
spectrum to the (nearly) zero expectation in the absence 
of a systematic error. To implement jackknife tests in 
the cross-spectrum framework, we calculate the cross- 
spectrum of differenced pairs of observations. More de- 
tails on the implementation can be found in LIO and 
Sll. We look for systematic effects in the following di- 
mensions. 

• Scan direction: We difference the data based 
on whether the telescope is scanning to the left 
or right. This test is sensitive to direction- or 
turnaround-dependent effects, such as detector 
time constants, or microphonics induced by tele- 
scope acceleration. 

• Azimuthal Range: We difference the data based 
on the azimuth at which it was observed. This is 
primarily a test for ground pickup. As discussed 
in Sll, we select the azimuth ranges for the dif- 
ference set based on the magnitude of observed 
ground emission on several degree scales. 

• Time: We have two time-dependent tests. First, 
we compare the first and second halves of the data 
for each field. Second, wc difference the bandpow- 
ers for 2008 from those of 2009. Unhke the other 
jackknife tests, this year- difference test is applied 
only to 150 and 220 GHz; there are no 95 GHz data 
from 2008. This year-difference is also implemented 
differently than the other jackknife tests. We calcu- 
late power spectra for each year and difference the 
bandpowers. Notably, this means that the test is 
less sensitive since the error bars must include sam- 
ple variance and beam uncertainties. These two 
tests are sensitive to long-timescale drifts in the 
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Fig. 1. — The results of the jackknife tests applied to the SPT 
data. Red squares denote the 95 GHz bandpowers, black circles 
mark the 150 GHz bandpowers, and the 220 GHz bandpowers are 
plotted with blue diamonds. Note that the 150 GHz uncertainties 
are smaller than the symbols. The top panel displays the results 
of the "first half - second half jackknife. The PTE to exceed 
the measured 0.90. The "left - right" jackknife is shown in 

the second panel and has a PTE of 0.21. The "azimuth-split" 
jackknife is plotted in the third panel and has a PTE of 0.84. 
Finally, the auto-spectra of the three frequencies are plotted in the 
bottom panel for comparison. 

calibration, pointing, detector properties, or any 
other time- varying effect. 

We find the results of each test to be consistent with 
zero. The results are plotted in Figure [T] Combining all 
three frequencies, the probabilities to exceed (PTE) for 
the "first half - second half , "left - right" and "azimuth- 
split" jackknives are 0.90, 0.21, and 0.84 respectively. 
The PTE of the observed x^'s is 0.50 for the three tests at 
95 GHz. For the four tests at 150 GHz and 220 GHz, the 
PTE are 0.48 and 0.61 respectively. There is no evidence 
for systematic effects in the SPT data. 

4. BANDPOWERS 

Applying the analysis described in j |2.3| to the 800 deg^ 
of sky observed with the SPT in 2008 and 2009 leads to 
the bandpowers shown in Figure [2] and tabulated in Ta- 
ble [l] The bandpowers, covariance matrix, and window 
functions are available for download on the SPT web- 
siteH 

A combination of the primary CMB anisotropy, fore- 
grounds, and secondary SZ anisotropics is detected at 
high S/N from £ = 2000 to 9400 in aU three frequency 
bands (see Figure [3| . The largest signal in the maps is 
the primary CMB anisotropy; a,t £ < 3000 the bandpow- 
ers trace out the Silk damping tail of the CMB power 
spectrum. The next largest signal is the CIB anisotropy 

^® http : / /pole. uchicago .edu /public/data/ reichardt 11/ 
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which becomes dominant around £ = 3000. At high mul- 
tipoles, the bandpowers follow the expected form for 
Poisson power from point sources. The spectral index 
of the point source contribution is a ~ 3.5 between 150 
and 220 GHz, characteristic of the Rayleigh- Jeans tail of 
a dust spectrum. Poisson fluctuations of radio sources 
become important as well at 95 GHz. Over these an- 
gular scales, the CMB and Poisson fluctuation power 
from point sources account for approximately 80% of the 
power in all frequency bands. We also detect the tSZ 
effect and clustering component of the GIB at high sig- 
nificance. 

5. COSMOLOGICAL MODELING 

We fit the SPT bandpowers to a model including 
lensed primary GMB anisotropy, secondary tSZ and kSZ 
anisotropics, and foregrounds. As stated in the previ- 
ous section, the GIB is the most important foreground 
and it is composed primarily of emission from a popula- 
tion of DSFGs. The power spectrum of these DSFGs can 
be divided into "Poisson" and "clustered" components. 
Radio sources are also important in the lower frequency 
bands. Finally, we include a template for the galactic cir- 
rus emission which is negligible by virtue of the selected 
sky coverage. 

Parameter constraints are calc ulated using the publicly 
available CosmoMCj^ package (Lewis & Bridle 2002). 
We have added two modules to nandle the high-i" data: 
one to model the foregrounds and secondary anisotropics 
and one to calculate the SPT and Planck/HFl likelihood 
functions. The code is a m odified and expand ed version 
of the modules discussed in Millea et al.| ( 2011 ) and used 
in Sll. These modules and instructions tor compiling 
them are available at the SPT website 

We include measurements of the primary GMB 
anisotropy, baryon acoustic oscillations (BAG), and Hub- 
ble constant (Hq) in all cosmological constraints pre- 
sented here. For the primary GMB, we inc lude both the 



seven -year WMAP data release (WMAP7, Larson et al. 
|2011|) and the recent large-scale SPT power spectrum 



We note that the same time-ordered data are used 
in this work and Kll. Hence we restrict the Kll data to 
non-overlapping multipoles {£ < 2000) . We also subtract 
a Poisson term from the Kll bandpowers to adjust for 
the different source m aski ng thresholds between the two 
works (for details see ^ 5.5 ). All parameter chains include 



low-redshift measurements of the Hubble constant H q us 
ing the Hubble Space Telescope (Riess et al. 2011) and 



the BAG feature using SDSS and 2dJ^'GRS data (Per 
cival et al.||2010D. Finally, we include Planck/HFlVTH 



measure ments at 217 GHz and in some cases 353 GHz 
( [Planck Gollaboration et al.||2011[ ). When comparing to 
the GIB model, we subtract the estimated Poisson radio 
contribution from the Planck/BFl GIB bandpowers. 

In some MGMG chains, we also compare the measured 
SZ power with external measurements of low redshift 
structure. We use local gal a xy clu ster abundances as 
measured by Vikhlinin et al. (2009'), which directly and 
tightly const rain cts and are consistent with o ther struc- 
ture probes dRozo et "aI1|2010| |Mantz et al.ir2 010 Van- 
' derlinde et al.p010| |Sehgal et al.p0il| . |Vikfilinm et airj 



■^^ http://cosmologist.info/cosmomc 

http: / /polo. uchicago.edu/public/data/rGichardtll/ 



( |2009| use the abundance of local (0.025 < z < 0.22) 
clusters to infer <Ts{n^/0.25 f-^'^ = 0.813 ±0.013 ±0.024, 
where the second set of errors is an estimate of the sys- 
tematic uncertainty due to the uncertainty in the masses 
of the clusters. This constraint will be referred to as the 
X-ray cts prior. 

5.1. Primary Cosmic Microwave Background 

Anisotropy 

We use the standard, six-parameter, spatially flat, 
lensed AGDM cosmological model to predict the primary 
GMB temperature anisotropy. The six parameters are 
the baryon density flth^, the density of cold dark matter 
ftch^, the optical depth to recombination r, the angular 
scale of the sound horizon at last scattering Gs, the am- 
plitude of the primordial density fluctuations ln[10^°As], 
and the scalar spectral index n^. 

We calculate the primary GMB power spectrum in the 
MGMG chains using PIGO (Fcndt & Wandch 20()7b]|a| 
instead of GAMB ( |Lewis et al. 1 ,2000) for speed. With 
lensing on, we find chain s run with PIGO (an d the fast 
WMAP likelihood code ( [Dvorkin fc h1I||2010[ ) modified 
for the latest WMAP7 likelihood) to be roughly 2000 
times faster than ones run with GAMB. PIGO allows us 
to replace the unlensed spectra plus "lensing template" 
approximation used in LIO and Sll with the full lens- 
ing treatment. We have trained PIGO using the January 
2011 version of GAMB with lensing including non-linear 
corrections turned on. Note that we observe no change to 
the foreground or secondary anisotropy constraints from 
the current data when using the former lensing template 
approximation. However, reducing the lensing power at 
£ = 3000 (by turning off lensing or neglecting non-linear 
effects) would increase the inferred kSZ power by an ap- 
proximately equal amount. For instance, neglecting non- 
linear effects would increase the kSZ power by ^ 1 /iK^. 

5.2. Thermal Sunyaev-ZeVdovich Anisotropy 

We adopt two different models for the tSZ power spec- 
trum. Sl l consi dered four tSZ models: a model by Seh- 



gal et al. (2010) which predicts high levels of tSZ power 
and three models which all predicted lower tSZ powers 



( Battagha et al.|2010[[Shaw et al.l2"oTo{|Trac et al.|2011[) 
Mere we have limited the model set to one from each 



group: the hi gh power model by Sehgal et al. (2010) and 
the model by Shaw et al. (2010) from the lower set. 



We use th e tSZ power spectrum pred icted by IShaw 
et~aL| (|20To| as the baseline model. , Shaw et al. ( |2010p 



investigate the impact of cluster astrophysics on the tSZ 
power spectrum using halo model calculations in com- 
bination with an analytic model for the IGM. These as- 
trophysical effects - in particular, the inclusion of non- 
thermal pressure support and A GN feedback - tend to 



suppress the tSZ p ower (see also Battaglia et al. 2011b 



Trac et al. 
that work 



20111). We use the baseline model from 
lereafter the Shaw model), which predicts 



Df^Q = 4.3 ^K^ at £ ^ 3000 and 152.9 GHz. The as- 
sumed cosmological parameters are {Qb, ^m, ^A, h, n^, 
o-g) = (0.044, 0.264, 0.736, 0.71, 0.96, 0.80). We use this 
model in all MGMG chains where another model is not 
explicitly specified. 

We also consider the tS Z power spectrum model pre- 
sented by ,Schgal et al. ( [20 10, hereafter the Sehgal 
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TABLE 1 
Bandpowers 









95 GHz 


150 GH 


z 


220 GHz 


i range 


4ff 


b {iiY?) a 




b (^K^) a 




b (mK2) 


o (mK2) 


2001 - 


2200 


2073 


211 A 


9.1 


212.9 


4.8 


285.2 


14.3 


2201 - 


2400 


2274 


138.4 


6.4 


130.1 


3.1 


203.0 


11.3 


2401 - 


2600 


2479 


107.5 


5.8 


97.4 


2.5 


188.3 


10.8 


2601 - 


2800 


2685 


75.1 


5.7 


73.9 


2.0 


171.7 


10.5 


2801 - 


3000 


2883 


51.9 


4.5 


55.1 


1.7 


175.6 


10.5 


3001 - 


3400 


3190 


47.5 


4.3 


41.5 


1.2 


171.2 


9.1 


3401 - 


3800 


3594 


32.8 


5.3 


32.5 


1.1 


179.8 


10.3 


3801 - 


4200 


3996 


35.5 


6.3 


32.5 


1.3 


208.2 


12.4 


4201 - 


4600 


4400 


37.2 


8.9 


35.9 


1.5 


224.0 


13.9 


4601 - 


5000 


4803 


14.8 


12.1 


37.5 


1.8 


273.3 


17.0 


5001 - 


5800 


5406 


31.3 


13.4 


45.4 


1.8 


341.8 


18.3 


5801 - 


6600 


6209 


38.6 


22.4 


57.0 


2.5 


422.0 


23.3 


6601 - 


7400 


7012 


26.8 


41.7 


70.6 


3.6 


455.1 


27.9 


7401 - 


8400 


7917 


-16.3 


78.3 


87.7 


5.0 


614.0 


39.8 


8401 - 


9400 


8924 


34.7 


156.0 


96.2 


7.3 


792.5 


62.1 




95 X 150 GHz 


95 X 220 GHz 


150 X 220 GHz 


2001 - 


2200 


2073 


207.6 


5.8 


207.2 


9.8 


226.0 


7.1 


2201 - 


2400 


2274 


129.6 


3.8 


131.7 


7.1 


148.6 


4.8 


2401 - 


2600 


2479 


94.0 


3.1 


97.7 


6.3 


112.8 


3.9 


2601 - 


2800 


2685 


72.1 


2.5 


80.7 


6.3 


97.4 


3.5 


2801 - 


3000 


2883 


49.9 


2.1 


59.5 


5.6 


78.6 


3.2 


3001 - 


3400 


3190 


37.9 


1.6 


39.3 


4.5 


69.2 


2.4 


3401 - 


3800 


3594 


22.9 


1.6 


22.4 


5.5 


64.8 


2.4 


3801 - 


4200 


3996 


25.0 


1.9 


38.3 


6.6 


72.4 


2.8 


4201 - 


4600 


4400 


25.5 


2.2 


23.3 


8.2 


77.6 


3.2 


4601 - 


5000 


4803 


24.1 


2.9 


42.0 


10.7 


94.5 


4.0 


5001 - 


5800 


5406 


31.4 


3.2 


53.6 


10.8 


111.4 


4.2 


5801 - 


6600 


6209 


38.5 


4.9 


51.9 


17.0 


139.3 


5.4 


6601 - 


7400 


7012 


48.1 


7.8 


40.0 


26.3 


169.9 


7.4 


7401 - 


8400 


7917 


57.0 


13.1 


122.5 


39.9 


211.2 


9.6 


8401 - 


9400 


8924 


83.9 


24.4 


109.3 


79.8 


282.0 


15.6 



Note. — £-band multipole range, weighted multipole value i?cff, bandpower _D, and bandpower uncertainty a for the six auto 
and cross-spectra of the 95 GHz, 150 GHz, and 220 GHz maps. The quoted uncertainties are based on the diagonal elements of 
the covariance matrix which includes the instrumental noise, Gaussian sample variance, and beam uncertainties. The sample 
variance of the tSZ effect and calibration uncertainty are not included. Calibration uncertainties are quoted in §2.2[ Point 
sources with flux > 6.4 mjy at 150 GHz have been masked out at all frequencies in this analysis. This flux cut substantially 
reduces the contribution of radio sources to the bandpowers, although DSFGs below this threshold contribute significantly to 
the bandpowers. 



model). Sehgal et al.| (2010) combined the se mi-analytic 
model for the intra-cluster medium (ICM) of Bode et al. 
( pC)09| ) with a cosmological N-body simulation to produce 
simulated tSZ and kSZ maps from which the template 
power spectra were measured. The assumed cosmolog- 
ical parameters are the same as with the Shaw model. 
At f = 3000, this model predicts Df^^ = 7 A fjK^ of tSZ 
power at 152.9 Ghz 



As described in SIX, the model of Shaw et al. (20101 
is also used to rescale both tSZ model templates as a 
function of cosmological parameters. Around the fiducial 
cosmological model, the tSZ amplitude roughly scales as 



h 
071 



1.7 

V0.80/ 



8.3 



0.044 



2.8 



(10) 



Both tSZ models exhibit a similar angular scale de- 



pendence over the range of multipoles to which SPT is 
sensitive. When the normalization of each model is al- 
lowed to vary, we detect similar tSZ p owe r for both (see 
Table [3]). However, as we discuss in ^7.2 the predicted 
model amplitude is crucial in interpretmg the detected 
tSZ power as a constraint on as and also leads to sig- 
nificant changes in the fit quality when the tSZ power is 
fixed to cosmologically-scaled model expectations. 

5.3. Kinetic Sunyaev-Zel'dovich Anisotropy 

Thomson scattering between CMB photons and clouds 
of free electrons that have a coherent bulk velocity pro- 
duce fluctuations in the observed brightness temperature 
of the CMB. This effect, known as the kinetic Sunyaev- 
Zel'dovich (kSZ) effect, leads to hot or cold spots, de- 
pending on whether the ionized gas is moving toward 
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Fig. 2.— The WMAP7 and SPT bandpowers. The SPT bandpowers at £ < 2000 are taken from Kll and are at 150 GHz only. At 
I > 2000, we show the bandpowers at 95, 150, and 220 GHz measured with the SPT in this work. Below I = 2000, the primary CMB 
anisotropy is dominant at all frequencies. On smaller scales, the GIB, radio sources, and secondary GMB anisotropics contribute to the 
signal. With the SPT source masking, the GIB is the largest source of power on sub-arcminute scales at 150 and 220 GHz. Due to the 
relative spectral behavior of the GIB and synchrotron emission, the 95 GHz bandpowers also have a significant contribution from radio 
sources. 



or away from the observer. The shift is identical to a 
change in the CMB blackbody temperature, except for 
tiny relativistic corrections. 

KSZ temperature anisotropy requires perturbations in 
the free electron density. This can be due to local pertur- 
bations in the baryon density, St, = Ph/Ph — ^ (where pi, 
is the mean baryon density in the universe), or ionization 
fraction, = Xe/xg — 1. Perturbations that are corre- 
lated with the large scale velocity field will produce an 
observable change in temperature of the CMB blackbody. 
The total contribution to the temperature anisotropy is 
given by the integral over conformal time. 



kSZ 



Tr 



CMB 



(ft) 



dr][a 



r{ri) 



Xe(77)]n-q, (11) 



where ctt is Thomson scattering cross-section, a is the 
scale factor, T{ri) is the optical depth from the observer 
to conformal time rj, n^fi is the mean electron density of 
the universe at the present day, and fi is the line of sight 
unit vector. Finally, 



(1 + (5,)(1 + Jb)^ 



(12) 



where v represents the peculiar velocity of free electrons 
at 77. 

The kSZ power spectrum can be broken down loosely 
into two components. The first is the kSZ signal from 
the post-reionization epoch, z < Zond, where we de- 
fine Zend as the redshift at which hydrogen reionization 
is complete. The second component comes from the 
epoch of reionization itself. Models of inhomogeneous 



(or patchy) reionization predict bubbles of free electrons 
around UV-emitting sources embedded in an otherwise 
neutral medium. Any large-scale bulk velocity of these 
bubbles will impart a temperature anisotropy onto the 
CMB. The kSZ signal from reionization is thus sensitive 
to fluctuations in both the gas density and ionization 
fraction. Note that in the case of instantaneous reion- 
ization, the post-reionization signal is the only source of 
kSZ power. 

We henceforth refer to the post-reionization compo- 
nent as the 'homogeneous kSZ' signal, and that from the 
reionization epoch as the 'patchy kSZ' signal. We now 
discuss our modeling of each in more detail. 

5.3.1. Homogeneous Kinetic Sunyaev-Zel'dovich Anisotropy 

For the homogeneous (or post-reionization) kSZ sig- 
nal, we adopt the non-radiative (NR) and cooling plu s 



star-fo rmation (CSF) models presented by Shaw et al. 



( |2011 |) 

anaiy 



These models are constructed by calibrating an 
ic model for the homogeneous kSZ power spectrum 
with two hydrodynamical simulations. The first simula- 
tion was run in the non-radiative regime, while the sec- 
ond included m etallicity-dep e ndent radiative cooling and 



star-formation. Shaw et al. (2011| measured the power 



spectrum of gas density fluctuations in both simulations 
over a range of redshifts, and used this to calculate the 
kSZ power spectrum for each case. They argue that these 
two cases are likely to bracket the true homogeneous kSZ 
power. For our fiducial cosmology (detailed in ^5.2) at 
£ = 3000, the NR template predicts 2.4 p.K'^, whereas the 
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Fig. 3. — The six auto- and cross-spectra measured with the 3-frequency SPT data. Overplotted on the bandpowers is the best-fit model 
for the fiducial set of model parameters. The bandpowers have not been corrected by the best-fit calibration or beam uncertainties in the 
MCMC chains; for reference, the best-fit temperature calibration factors at 95, 150, and 220 GHz are 0.999, 0.997, and 1.003 respectively. 
In addition to the complete model (black lines), each individual model component is shown. The tSZ effect is marked with the blue solid 
line. The best-fit kSZ power is near-zero and off-scale. The Poisson power from DSFGs and radio galaxies are shown by solid orange 
and green lines respectively. The clustered component to the DSFGs is shown with a orange dot-dash line. 



CSF simulation predicts 1.6 jiK^ . The scaling of these 
models with cosmological parameters is given approxi- 
mately by 




Note that, for both models, we assume that reioniza- 
tion occurs instantaneously at Zend = 8. Both models 
also assume that Helium remains neutral; if Helium is 
singly ionized at the same time as hydrogen, the pre- 
dicted power would increase by 16%. 

5.3.2. Kinetic Sunyaev-Zel'dovich Anisotropy from Patchy 
Reionization Scenarios 

For the patchy kSZ signal, we adopt the model pre- 
sented by Zll and briefly summarized here. This model 
starts from a matter simulation. For every position, the 
linear matter over-density is calculated within some ra- 
dius. The halo collapse fraction is estimated from this 
over-density, and translated into an expected number of 
ionizing photons based on the number of collapsed halos 
above the atomic cooling threshold mass. The number 



of ionizing photons is then compared to the number of 
hydrogen atoms within this sphere. If there are sufhcient 
ionizing photons, the sphere is labeled ionized. If not, a 
smaller radius is set and the algorithm repeated until the 
resolution of the simulation box is reached. The resulting 
ionization field is used in conjunction with the underly- 
ing density and velocity fields to compute the patchy kSZ 
power spectrum. 

As discussed in Zll, the detailed reionization history 
used for the kSZ template is unimportant, since the tem- 
plate shape is robust to the duration and the mean red- 
shift of reionization. The amplitude of the patchy kSZ 
signal is nearly proportional to the duration of reion- 
ization (with a mild rcdshift dependence), making it an 
excellent probe of reionization. The kSZ effect is comple- 
mentary to the large scale CMB polarization anisotropy 
(the reionization bump) that provides a constraint on 
the global timing of reionization, and to QSO and Lya 
emitter observations that study the very end of reioniza- 
tion. For the fiducial model template provided by Zll 
and used in this work, reionization begins at z = 11 and 
concludes at z = 8. 

5.4. Cosmic Infrared Background 

The CIB is produced by thermal emissi on from DSFGs 
over a very broad range in redshift ( |Lagache et al.|[2005 
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Marsden et al. 20091. The dust grains, ranging in size 
irom a tew molecules to 0.1 mm, absorb light at wave- 
lengths smaller than their size, and re-radiate it at longer 
wavelengths. The result of this absorption is roughly 
equal amounts of energy in the CIB as in the unpro- 
cessed s tarlight that m akes up the optical/UV back- 
ground pwek fc Aremlt 1998; Fixsen et al.||1998[ ). 

The frequ ency dependence of the CIB (and radio galax- 
ies in [5.5) is naturally discussed and modeled in units 
of flux density ( Jy) rather than CMB temperature units. 
However, as the power spectrum is calibrated in CMB 
temperature units, we need to determine the ratio of 
powers between different frequency bands in CMB tem- 
perature units. To calculate the ratio of power in the 
Vi X Vj cross-spectrum to the power at (the arbitrary fre- 
quency) uq in units of CMB temperature squared, we 
multiply the ratio of flux densities by 



dB I dB I 
dT I '^0 dT 1 



dB I 
dT I 



dB I 
: dT I 



(14) 



where B is the CMB blackbody specific intensity evalu- 
ated at TcMB, and I'i and I'j are the effective frequencies 
of the SPT bands. Note that i^i may equal i^j. The ef- 
fective frequencies of the SPT bandpasses for each fore- 
ground are presented in {5.7 



5.4.1. Poisson term 

Statistical fluctuations in the number of DSFGs lead to 
a Poisson distribution on the sky. DSFGs are effectively 
point sources at these angular scales, so they contribute 
a constant Ci (thus Dg oc £^) to the SPT bandpowers. 
The bulk of the millimeter-wavelength CIB anisotropy 
power is expected to come from faint (< 1 mJy) DSFGs 
- so the residual CIB power is largely independent of the 
source mask flux threshold. This is a key assumption 
in combining the Planck/HFl and SPT data as the flux 
cuts differ by a factor of eight at 220 GHz. The 217 GHz 
Planck/RFl cut is at 160 mJy while the 150 GHz SPT 
flux cut of 6.4 mJy scales to approximately 22mJy at 
220 GHz. We test the im pact of the flux c u t with the 
modeled source coun ts of Bethermin et al. (2011) and 
Negrello et al. ( 2007 1 , and find the two models predict a 
difference in the Poisson power of 4% or 16% respectively. 
Hence in the model fitting, we scale the modeled DSFG 
Poisson power up by 10% when comparing the model to 
the Planck/HFl CIB bandpowers. We have tested chang- 
ing this assumption by tens of percent and find the cos- 
mological constraints are insensitive to this - the Poisson 
term is subdominant at all multipoles {£ G [80,2240]) in 
the Planck/KFl bandpowers. 

The Poisson power spectrum of the DSFGs can be writ- 
ten as 



^t,vi,v2 ~ -^3000 ^i^l: ''2 



i 



3000 



(15) 



where ij^, encodes how t he CIB brightness scales with 
observing frequency (see ^ 5.4.3). ffooo ^^'^ amplitude 
of the Poisson DSFG power spectrum at ^ = 3000 and 
frequency i/q. Note that we neglect the decohering effect 



of frequency dep endence varying from source to source 
Hall et al. (2010') argued that the spectral index variation 
in the DSFG population was less than aa = 0.7. This 



level of variation has a negligible effect on the inferred 
Poisson powers across the SPT frequency bands. 

5.4.2. Clustered term 

DSFGs (and galaxies in general) trace the mass distri- 
bution, so they are spatially correlated. This leads to a 
"clustering" term in the DSFG power spectrum in addi- 
tion to the Poisson term. On the physical scales of rele- 
vance, galaxy correlation functions are surprisingly well- 
approximated by a power law. The clustering shape we 
consider is thus a phenomenologically-motivated power 
law with index chosen to match the correlation pro per- 
ties observed for Lyman-break galaxies at z ~ 3 (|Gi 



avalisco et al.| |1998] [Scott fc White] |1999[ ). This form 



as been used previously by a number of millimeter 



and submillim eter-waveleng t h experiments ( Viero et al 



(12009 ) ; Dunklc y et al. jgOll ) ^lanck Collaboration efaT 
pbll| ); lAddison et aI.| " pDlT| , Sll). We adopt this 
power law" model to represent the non-linear regime, 
and use a clustering template of the form Z?^ c>c £'~''^ in 
all ca ses u nless specified. 

In §9.1[ we consider two alternate forms for the clus- 
tering term. First, we free the power-law exponent, i.e., 
I?^ oc with 7 free . Seco nd, we add the linear-theory 
model by Hall et al. (20101 with a free amplitude to the 
power-law clustering term. The fiducial model in that 
work was calculated assuming that light is a biased tracer 
of mass fluctuations calculated in linear perturbation the- 
ory, with redshift and scale-independent bias. Further, 
the model assumed all galaxies have the same greybody 
SED, and the redshift distribution of the luminosity den- 
sity was given by a parametrized form, with parameters 
adjusted to fit SPT, BLAST, and Spitzer power spec- 
trum measurements. We find that the angular template 
shape depends only slightly on the observing frequency 
and choose to use the 150 GHz model spectrum. Here we 
refer to this power spectrum shape as the "linear theory" 
template. As will be discussed in S|9j the data strongly 
prefer the power-law form over the linear theory clus- 
tering template. The data are also consistent with the 
assumed power law exponent of 7 = 0.8. 

5.4.3. Frequency dependence 

Unless otherwise noted, all results presented in this 
work assume the modified black-body (modified BB) 
model for the frequency dependence of the CIB: 



77, = iy^B,{T), 



(16) 



where Bi, (T) is the black-body spectrum for temperature 
T, and /3 is an effective dust emissivity index. T and /3 
are free parameters with uniform priors T e [5, 35 K] 
and /3 e [0, 2]. The upper limit on the temperature prior 
is unimportant because the SPT frequencies are in the 
Rayleigh- Jeans tail of >15 K objects. Therefore, there 
would be a perfect degeneracy between 30 K and 100 K 
dust. We have tested instead fixing the t emperature to 
9.7 K with a wider (3 range as was done in [Addison et al. 
(2011) and found little effect on the results reported in 
this work. 

For some results , we also exami ne the single-SED 
model presented by Hall et al. (2010) (hereafter the Sin- 
gle SED model); however, we do liot find substantive 
differences in the results from the modified BB model. 
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This model assumes a single SED for all DSFGs. This 
SED is a modified black-body with T = 35 K and /3 = 2. 
The DSFGs have a density as a function of redshift 
parametrized by a Gaussian with mean Zc and variance 
cr^. The number counts are set to zero at z < 0. We 
fix all but one par ameter in th e mod el; we allow to 
vary freely. As in Hall et al. (2010), we set ctz = 2. 



The frequency dependence predicted by the Single SED 
model has a slight ^-dependence; we ignore this and use 
the frequency scaling at £ — 3000. We use the model at 
£ = 3000; it has a slight ^-dependence which we ignore. 
The Single SED model has one fewer parameter and is 
more easily interpretable physically than the modified 
BB model. However, the data are in mild tension with 
the model - the preferred Zc is near the prior cutoff at 0. 
This does not degrade the quality of the fit to the data, 
but it leads to a one-sided distribution of GIB spectral in- 
dices which results in a slight artificial reduction in some 
parameter uncertainties. 

A simplifying assumption we make for the analyses 
presented here is that the various components of the 
GIB (Poisson and clustering terms) have the same fre- 
quency dependence. Even in the context of the Single 
SED model, we expect their spectral indices to differ to 
some degree, because the weighting of sources as a func- 
tion of redshift will be different for the two terms. Our 
assumption here is, effectively, that the redshift distribu- 
tions of the source populations giving rise to the Poisson 
and clustered components are not sufficiently different as 
to cause a significant difference in frequency dependence. 
We have tested allowing separate frequency dependencies 
for the Poisson and clustered GIB terms and find mini- 
mal impact on the SZ constraints. 

We have also checked whether the higher fre quency, 
217 and 353 GHz Planck/R Fl GIB bandpowers ( ,Planck 
Gollaboration et al. 2011) are consistent with the as- 
sumption ot a single frequency dependence for the Pois- 
son and clustered terms. If the frequency dependence 
of each term was different, then we would expect to see 
variations in the GIB power ratio between 217 GHz and 
353 GHz as a function of angular scale since the rela- 
tive power in each term is ^-dependent. We find that the 
GIB power ratio is consistent with a single spectral in- 
dex with a2i7/353 = 3.22 ±0.08. With (9-1) dof, the nuU 
hypothesis of a single spectral index has = 8.0 - the 
probability to exceed this is 0.44. This spectral in- 
dex is flatter than preferred by the lower frequency SPT 
data; however, we expect the GIB frequency dependence 
to flatten towards higher frequencies as the rest frame 
emission frequency approaches the peak of the grey body 
dust emission. 

5.4.4. tSZ-CIB correlations 

As tracers of the same dark matter distribution, we 
expect some degree of correlation between the tSZ and 
GIB power spectra. In previous papers, we have ar- 
gued that the correlation is likely to be small since most 
of the DSFGs contributing to the GIB are in the field 
rather than in the massive galaxy clusters that contribute 
most of the tSZ power. Unfortunately, observational con- 
straints on the IR flux (and therefore tSZ-GIB correla- 
tion) are largely limited to low-redshift, high-mass clus- 
ters. The observations are also at frequencies close to the 



peak of the GIB emission, introducing a significant scal- 
ing uncertainty in translating the results to the longer 
wavelengt hs irnportan t for tSZ surveys. The simulations 
of Sehgal et al. ( |2010 1 show that a model with significant 
correlations at 148 GHz can be constructed without vi- 
olating current observational constraints by concentrat- 
ing the GIB galaxies in the lower-mass, higher-redshift 
galaxy clusters. 

We adopt a simplified treatment of tSZ-GIB correla- 
tions in this work. We assume the correlation, ^, is inde- 
pendent of angular multipole, and that both the Poisson 
and clustering GIB components correlate with tSZ fluc- 
tuations. Note that the choice of which GIB terms corre- 
late does not significantly change the resulting tSZ and 
kSZ constraints, although it does rescale the amplitude 
of the inferred correlation. 

The tSZ-GIB correlation introduces a term to the 
model of the form: 



D 



tSZ-CIB 

' -Ml ^^2 



tsz 



1-^1,1^2-1^2) V \-L^e,vo,u2-^e,i^i-i^i- 

(17) 

Here Df;"-;:^ is the power due to correlations (which 
should be negative below the tSZ null of 217 GHz, i.e. 
C < 0), Df^^^^^ is the tSZ power spectrum, and 
is the sum of the Poisson and clustered GIB components. 

5.5. Radio galaxies 

High-flux point sources in the SPT maps are coincident 
with known r adio sources from lo ng- wavelength catalogs 
(e.g. SUMSS, p"auch et al.|r2003[ ), and hav e spectral in- 
dices consisteiit with synchrotron emission ( Vieira et al. 
20101). The Poisson radio source term can be expressed 



as: 



„r ar+0.51n())„j^,„2)'''7- 



3000 



, (18) 



where rj^ 

of the spectrum to the base frequency. i^Jooo 



is the ratio of the frequencies 
is the am- 
plitude of the Poisson radio source power spectrum at 
£ — 3000 and frequency vq. The mean spectral index of 
the radio source population is with scatter ar- 

The radio source power (£^3000) map depends 

approximately linearly on the flux above which discrete 
sources are masked. This reflects the fact that S'^dN/dS 
for synchrotron sourc es is nearly independent of S (e.g. 
de Zotti et al. ( |2005 )). As a result, the Poisson radio 
source power m the Kl 1 bandpowers (which had a higher 
flux cutoff) is expected to be 9.2 ± 2.4 fiK^ higher than 
in the 150 GHz bandpowers presented here. We correct 
for this by subtracting a Poisson term with amplitude 
Z?3gQQ = 9.2 /iK'^ from the Kll bandpowers in the pa- 
rameter fits. Note that the restriction to £ < 2000 we 
impose on the Kll bandpowers renders the uncertainty 
on the Poisson power insignificant. We have tested vary- 
ing the subtracted Poisson power and found no impact 
on the results. 

Following the treatment in Sll, we apply a strong prior 
to the Poisson radio source term due to sources below the 
SPT detection threshold. First, we fix the mean spectral 
index (ar = —0.53) and sc atter (ct,. — 0). Second, based 



on the de Zotti et al. ( 2005 ) source count model, we apply 
a Gaussian prior to the 150 GHz amplitude of 1.28 ± 
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0.19 fiK^. We have tested the impact of removing the 
radio priors. The data prefer parameter values consistent 
with the priors, and there are no significant changes to 
the CIB and SZ constraints. The strong radio priors are 
included in all cas es presented here. 



As discussed in [Hall et al. (20101, the clustering of 
radio galaxies at Sr'l' frequencies is likely to be negligi- 
ble. Clustering on scales of a few arcminutes is a ~ 5% 
modulation of the mean, and the intensity of the radio 
source mean is quite low. Furthermore in the unified 
model of active galactic nuclei, the dominant sources at 
> 90 GHz (blazars and flat spectrum radio quasars) are 
active galactic nuclei with their jets aligned towards the 
observer. The effect of random alignment would serve 
to further dilute the clustering signal. Supporting this 
argument are conservative extrapolati ons of upper limit s 
on arcminute scale power at 30 GHz (Sharp ct al. 2010). 
If entirely due to radio point source power, the 30 GHz 
limit translates to an upper limit of 0.1 fiK^ at 150 GHz. 
We do not include a clustered radio component in any 
results in this work. 



5.6. Galactic Cirrus 

Galactic cirrus is the final (and smallest) foreground 
we incl ude in our modeling. Following the cirrus treat- 
ment in Hall et al.| ([2010) and Kll, we cross-correlate the 
SPT maps with the mod el 8 galactic dust predictions of 



Finkbeiner et al.] (1999). Galactic cirrus is detected in 
all frequency bands at a significance of ^ 3cr with an 
angular dependence of: 



T^CIV _ 7-)Cir,!^i,l/2 

^£,1/1,1^2 ~ ^3000 



3000 



-1.2 



(19) 



The measured powers are 0.16, 0.21, and 2.19 ^K^ at 95, 
150, and 220 GHz. Note that the 150 GHz amplitude is 
different from the value quoted in Kll because of the 
different field weights in the two analyses. We fix the 
amplitude of the galactic cirrus in each frequency band; 
allowing it to vary with a prior based on the measurement 
uncertainty has no impact on any parameter constraints. 

5.7. Effective frequencies of the SPT bands 

Throughout this work, we refer to the SPT frequen- 
cies as 95, 150, and 220 GHz. The data are calibrated to 
CMB temperature units, hence the effective frequency is 
irrelevant for a CMB-like source. However, the effective 
band center will depend (weakly) on the source spec- 
trum for other sources. The actual SPT spectral band- 
passes are measured by Fourier transform spectroscopy, 
and the effective band center frequency is calculated for 
each source, assuming a beam-filling source with a nomi- 
nal frequency dependence. The bandpasses differ slightly 
from 2008 to 2009 due to the focal plane refurbishment. 
We calculate the effective frequency for each year's spec- 
tral bandpass and ave rage them according to the yearly 



weights presented in S2.1 For an a = —0.5 (radio-like) 
source spectrum, we find band centers of 95.3, 150.2 and 
214.1GHz. For an a = 3.5 (dust-like) source spectrum, 
we find band centers of 97.9, 153.8 and 219.6 GHz. The 
band centers in both cases drop by ~ 0.2 GHz if we reduce 
a by 0.5. We use these radio-like and dust-like band cen- 
ters to calculate the frequency scaling between the SPT 



bands for the radio and DSFG terms. For a tSZ spec- 
trum, we find band centers of 97.6, 152.9 and 218.1 GHz. 
The ratio of tSZ power in the 95 GHz band to that in 
the 150 GHz band is 2.78. The 220 GHz band center is 
effectively at the tSZ null; the ratio of tSZ power in the 
220 GHz band to that in the 150 GHz band is 0.00004 for 
a non-relativistic tSZ spectrum. 

6. RESULTS 

6.1. Baseline Model Results 

We adopt a baseline model with the six ACDM pa- 
rameters and six additional parameters for the secondary 
CMB anisotropics and foregrounds described in |5] Four 
of these extra parameters represent the power in the tSZ 
effect, kSZ effect, DSFG Poisson term, and DSFG clus- 
tering term. The final two parameters describe the fre- 
quency dependence of the DSFG anisotropy. All DSFG 
terms are assumed to have the same frequency scaling, 
based on the modified BB model. As templates for the 
angular shape of tS Z and kSZ power spectra, we u se the 



Shaw et al. (2010) tSZ model and Shaw et al. (2011) 



Ctih' kSZ mbdeTThe DSFG Poisson power (like aTHnT 
resolved Poisson power) has the shape Di cx £^ while 
the DSFG clustering term is descr ibed b v the power-law 
De cx template presented in { 5.4.2 Figure Is] plots 
the best-fit values for each component on top of the SPT 
bandpowers. 

Table[2]shows the improvement in the quality of the fits 
with the sequential introduction of free parameters to the 
original ACDM primary CMB model. Adding additional 
free parameters beyond the tSZ amplitude does not sig- 
nificantly improve the quality of the fits. However, small 
improvements are seen with the introduction of tSZ-CIB 
correlations or treating the clustered and Poisson fre- 
quency scalings as independent. 

The additional parameters discussed in ^Js] are either 
fixed or constrained by strong priors. A Poisson popu- 
lation of radio galaxies is included with a fixed spectral 
index ar = —0.53, zero intrinsic scatter, and an ampli- 
tude of £'3QQQ = 1.28 ± 0.19 /iK^. Weakening these radio 
priors - removing the amplitude prior and setting a uni- 
form prior on the radio spectral index of G [^1-5,0] 
- does not change the inferred SZ levels. We always 
include the constant Galactic cirrus power described in 
|5.6[ We assume no correlation between the tSZ effect 
and CIB in the baseline model. This assumption is re- 
laxed in j |6.2.lj The for the 90 SPT bandpowers (84 
dof since the ACDM parameters are essentially fixed by 
external data and there are six model parameters be- 
yond ACDM) is 73.3; the PTE for this is 79%. We 
have tested fixing the ACDM parameters parameters to 
the best-fit values, and find the derived constraints on 
the six extended model parameters are unchanged. This 
baseline model fits the SPT data well and provides the 
most concise interpretation of the data. 

SZ parameter constraints for the baseline model and 
extensions are listed in Table [3] and are presented in 
Figure H For the tSZ power, we find D^g^o = 3.65 ± 
0.69 /iK^ which is consistent with the predictions of the 
Shaw model. We find a strong 95% CL upper-limit on 
kSZ power at -D|ooo < 2.8 /LtK^, although we stress that 
this limit is dependent on our DSFG modeling assump- 
tions. We have applied a positivity prior to the kSZ 
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TABLE 2 
Delta for model components 



term 


dof 


Ax' 


CMB 


6 


(reference) 


DSFG Poisson 


3 


-1645 


DSFGs NL Clustering- 


1 


-218 


Radio Poisson 





-164 


tsz 


1 


-55 


kSZ 


1 





Single SED 


-1 





Ctp / 


2 


-4 


tSZ-CIB correlation 


1 


-2 


P clustered DSFGs 


1 





LT clustered DSFGs 


1 






Note. — Improvement to the best-fit x ^ additional 
terms are added to the model. Terms above the double line 
are included in the baseline model, with each row showing the 
improvement in likelihood relative to the row above it. For 
rows below the double line, the Ax' is shown relative to the 
baseline model rather than the row above it. Here we have 
assumed the modified BB frequency scaling of the GIB except 
for the row labeled "Single SED" which assumes the alternate 
Single SED frequency scaling. The row labeled "tip 7^ Qc" 
allows the Poisson and clustered power from DSFGs to have 
independent frequency scalings. The row labeled "tSZ-GIB 
correlation" introduces a correlation coefficient between the 
tsz and GIB power spectra. The last two rows add freedom 
to the clustered DSFG model shape, either by varying the 
power-law exponent ( clustered DSFGs" ) or adding a free 
amplitude linear theory template ("LT clustered DSFGs"). 



power; without this prior, the median kSZ power is neg- 
ative by 1.2(7. The SZ constraints are independent of 
the assumed template shape; we find similar power con- 
straints for all tSZ and kSZ model templates considered. 

The current tSZ and kSZ constraints are consistent 
with earlier measurements of the SZ power. Sll mea- 
sured D^g^o + O.bD^^^Q = 4.5 ± 1.0 /iK^. The equivalent 
constraint from these data is {1.03A)Dlf^Q + O.5D3000 = 
4.27 ± 0.58 /iK^. The multiplicative factor in parenthe- 
ses before the tSZ term corrects for the different effective 
bandcenters in the two datasets 



Dunkley et al. (20111 



measured -Dsooo + ^aooo = 6.8±2.9 /LtK^, while these data 
prefer {l.l7)Dll^o + ^3^00 = 5.27 ± 0.76 mK^ We note 



that the two analyses differ in that Dunkley et al. (2011 1 
fixed the ratio of thermal and kinetic SZ power, while 
this ratio is free in this work. The three datasets show 
excellent consistency. 

We detect both the DSFG Poisson and clustering 
power with high significance, with powers of i'fooo ^ 
7.54 ± 0.38 /iK2 and i:»|ooo = 6.25 ± 0.52 ^K^ respec- 
tively at 150 GHz. The effective spectral index from 
150 to 220 GHz of the DSFG power is found to be 
a = 3.56 ± 0.07. The DSFG constraints are discussed 
in more detail in ^J9l and are in line with both theoreti- 
cal expecta tions and previous work (HIO; Sll; Dunkley 
eFaI][20lT| ). 



6.2. Results for Alternative Models 

Here we explore modifications to the baseline model 
and the impact on SZ and CIB constraints. First, we 



consider allowing a correlation between the tSZ effect 
and CIB, and show the impact on the derived SZ power. 
The SZ interpretation is robust to the other CIB model 
variants we examined. Second, we examine model ex- 
tensions to the ACDM cosmology, and show they do not 
affect the SZ and CIB parameters. 

6.2.1. tSZ-CIB Correlations 

We first consider the effect of a correlation between 
the tSZ and CIB as prescribed by equation [17] Note 
that DSFG over-densities in galaxy clusters would lead 
to an anti-correlation since the tSZ effect causes a flux 
decrement at the SPT frequencies. The correlated power 
is highly degenerate with the kSZ and tSZ power as 
shown in Figure [5] We see some dependence on the 
kSZ template shape (not shown); the patchy template, 
which is falling instead of rising slightly across the an- 
gular multipoles of interest, is less degenerate with the 
correlation coefficient than the CSF homogeneous kSZ 
model. Increasing anti-correlation increases the allowed 
kSZ power, while reducing the tSZ power, until ^ ~ —0.4. 
The slope of the trade-off between ^ and kSZ power ap- 
proaches zero for very negative values of ^. This is be- 
cause the correlated power is proportional to the tSZ 
power, and the tSZ power approaches zero in this limit. 

Decreasing tSZ power with increasing anti-correlation 
contradicts our naive expectations. In single-frequency 
bandpowers (e.g. only 150 GHz), increasing tSZ-CIB 
anti-correlations would increase the allowed tSZ contri- 
bution for any observed power level. However, this pic- 
ture changes when we consider the other SPT frequency 
combinations. A tSZ-CIB anti-correlation reduces the 
power at 150 x 150 (which could be compensated by in- 
creasing tSZ power) but it also has ^50% larger effect 
at 150 X 220 (where there is effectively zero tSZ power). 
Conversely, a tSZ-CIB anti-correlation reduces the power 
at 95 X 150 by a similar amount to 150 x 150, whereas 
there is more tSZ power in the 95 GHz band. Adding 
kSZ power more effectively cancels the tSZ-CIB correla- 
tion term than adding tSZ power in the most sensitive 
combinations of the three frequencies. 

Due to this degeneracy, introducing a tSZ-CIB corre- 
lation substantially degrades the constraints on both the 
tSZ and kSZ power. In a model with ^ as a free parame- 
ter, the measured tSZ power is D^^^ = 3.26 ± 1.06 /^K^. 
The 95% confidence upper limit on the kSZ power is now 
^3000 < 6-7 /iK^. The SPT data better constrain the hn- 
ear combination of Df^^ + ^-^D^mo = 4-60 ± 0.63 //K^, 
although the limiting case with zero tSZ and high kSZ 
power is disfavored. The expanded parameter space re- 
sembles the constraint in previous SPT work (LIO, Sll); 
the additional frequency information in this analysis is 
constraining the CIB model instead of breaking the de- 
generacy between tSZ and kSZ power. To facilitate com- 
parison we express the constraint in the same basis as 
Sll; we find a constraint of (1.034)^*32^ -f ^.hD'^l^ = 
4.71 ± 0.64 /iK^. The multiplicative factor in paren- 
theses before the tSZ term corrects for the different 
effective bandcenters in the two datasets. Sll found 
^3000 + 0.5L»^^o^o = 4.5 ± l.O/xK^. The two measure- 
ments are completely consistent with each other. We 
will discuss the inferred correlation and its impact on 
CIB constraints in ^ 
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TABLE 3 

SZ CONSTRAINTS 



tSZ Model 


kSZ Model 


CIB model 


^fooo (mK') 




Shaw 


CSF 


modified BB 


3.65 ±0.69 


< 2.8 


Shaw 


CSF 


modified BB w. correlations 


3.26 ± 1.06 


< 6.7 


Shaw 


patchy 


modified BB w. correlations 


3.66 ±0.86 


< 5.7 


Shaw 


fixed CSF 


modified BB 


3.34 ±0.58 


(1.57) 



Note. — Measured tSZ power and 95% confidence upper limits on the kSZ power at ^ = 3000 in the SPT 150 GHz band. 
Two tSZ model templates (the Shaw and Sehgal models) and two kSZ models (the CSF homogeneous kSZ model and the patchy 
reionization kSZ model) have been considered. In most cases, we show results only for the Shaw tSZ template and CSF kSZ 
template as the other cases are essentially identical - indicating the data are insensitive to the differences in the modeled angular 
dependencies. In all cases, the modified BB model is assumed for the CIB frequency scaling. As an extension to the CIB model, 
we present constraints when correlations are allowed between the tSZ and CIB. This correlation allows twice the kSZ power and 
weakens the tSZ power constraint by 50%. The results depend somewhat on the kSZ template when tSZ-CIB correlations are 
introduced. The real kSZ spectrum should be an mixture of the two kSZ templates. Finally, we show the tSZ constraint derived 
when fixing the kSZ power to that predict ed by the CSF homogeneous kSZ model. The CSF kSZ model predicts 1.57/^K^ for 
the WMAP7 cosmology detailed in §5.3. 1| and is scaled with cosmology according to the prescription in that section. 
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Fig. 4. — 2D likelihood surface for the tSZ and kSZ power in the SPT 150 GHz band at £ = 3000. 1, 2, and 3-0- constraints are 
shown in red, yellow, and blue respectively. The left panel shows the constraints for the modified BB CIB frequency scaling. The right 
panel displays the likelihood surfaces for the modified BB CIB frequency scaling with the introduction of free parameter, ^, describing the 
correlations between the tSZ and CIB power spectra. The tSZ-CIB correlation is degenerate with the kSZ effect, thereby weakening the 
kSZ limits. 



6.2.2. Beyond ACDM 

Lastly, we consider whether the derived SZ and CIB 
powers would change with an alternate cosmological 
model. Current measurements of the primary CMB 
anisotropy contain hints of parameters beyond the six 
ACDM parameters (Rcichardtet al.|2009 Komatsu et al. 
2011||Dunkley at al.||2011( Kll). T'he evidence tor these 
model extensions is not yet compelling; however, it is in- 
teresting to test whether the extended models affect the 
SZ constraints. We consider three extensions: running of 
the scalar spectral index of primordial fluctuations, mas- 
sive neutrinos, and extra neutrino (or other relativistic 
particle) species. We find no evidence for a degeneracy 
between any extended parameters and the kSZ or tSZ 
power. The kSZ and tSZ constraints are unchanged un- 



der these three ACDM model extensions. 

We also consider using an unlensed primary CMB 
power spectrum. We see a significant increase in the 
kSZ power when lensing is turned off. The inferred kSZ 
power increases by about the same amount as the pri- 
mary CMB decreases in the absence of gravitational lens- 
ing at £ ^ 3000. This produces a significant, ^ 
shift in the kSZ power. However, current CMB data rule 
out no-lensing scenarios at more than 5cr (Kll). The 
induced uncertainty on the kSZ level for the current con- 
straint on the lensing amplitude, A^^^^ = 0.95 ± 0.15, is 
only 0.4 /iK^. This is much smaller than the Icr uncer- 
tainty on kSZ, and therefore can be neglected. 



7. THERMAL SZ INTERPRETATION 
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Fig. 5. — 2D Likelihood curves for the correlation between the tSZ and CIB spectra versus the kSZ power {Left panel) or tSZ power 
{Right panel) in fiK^. The filled contours show the 1, 2, and 3 a constraints. The kSZ template consists of the sum of CSF homogeneous 
kSZ model and patchy kSZ contribution for a model where reionization begins at 2 = 11 and ends at z = 8. Increasing anti-correlation 
increases the allowed kSZ power and reduces the tSZ power. The data prefer anti-correlation suggesting that DSFGs are over-dense in 
galaxy clusters. 

We now focus on comparing our measurements of the 
tSZ power spectrum with the predictions of recent mod- 
els and simulations. We combine measurements of the 
primary CMB power spectrum, Hq, and BAO with the 
tSZ signal to measure as and the sum of the neutrino 
masses, m^. 

7.1. Thermal SZ constraints 

The models for the tSZ power spectrum are based on 
simulations with a fixed cosmology. Therefore, to cor- 
rectly compare the measured tSZ power with theoret- 
ical predictions, we must rescale the theoretical power 
spectra so that they are consistent with the cosmological 
parameters at each point in the MCMC chains. Recall 
that the cosmological parameters are derived solely from 
constraints on BAO, Hq, and the primary CMB power 
spectrum and so are not directly influenced by the am- 
plitude of the measured tSZ power. 

Following Sll, we define a dimensionless tSZ scaling 
parameter, Atsz, where 



ntsz 

^3000 



(20) 



DgpQQ is the measured tSZ power at ^ = 3000 and 
$3QQg(Ki) is the template power spectrum normalized to 
be consistent with the set of cosmological parameters '/Ci' 
at step i in the chain. We obtain the preferred value 
of Atsz by combining the measurement of I?3ooo with 
the cosmological parameter constraints provided by the 
primary CMB measurements. As described in Section 
|5.2[ we consider two template models for the tSZ power 
spectrum, the 'Shaw' and 'Sehgal' models. These two 
templates bracket the range of tSZ models considered by 
Sll. 



To determine the cosmological scaling of $3QQQ(Ki), we 
run the Shaw model at each step in the chain. The 
approximate cosmological scaling of the Shaw model 
around the fiduci al m odel cosmological parameters is 
given by Equation |10[ 

In the right panel of Figure [6j we show the constraints 
obtained on Atsz • The solid black line in each plot shows 
the results for the baseline model using the Shaw tSZ 
template, the CSF kSZ template, and the modified BB 
CIB frequency scaling, while the green line shows those 
obtained when a free tSZ-CIB correlation is allowed in 
our modeling. The blue dashed line shows the results ob- 
tained using the Sehgal tSZ model instead. The numer- 
ical constraints are also given in Table |4j Recall that a 
value of Atsz — 1 implies that the measured tSZ power is 
consistent with the tSZ template predictions, while Atsz 
less (greater) than 1 indicates that the model overesti- 
mates (underestimates) the amplitude. Note that sample 
variance will cause deviations from unity on the order of 
0.04. 

The preferred value of Atsz for the Shaw template is 
1.2 (T below unity (from the constraints given on In Atsz), 
while the preferred value of Atsz for the Sehgal template 
is 2.9(7 below unity. Although both templates predict 
more tSZ power than is measured, the Sehgal template 
is in significant tension with the data. As an aside, we 
note that using the Single SED CIB frequency scaling 
instead of the modified BB has little effect on the tSZ 
amplitude. Including a tSZ-CIB spatial correlation shifts 
the preferred value of Atsz to a slightly lower value, but 
principally increases the error bars (by roughly 33%). 

In Table [5] we evaluate the goodness of fit of our SZ 
models by running new MCMC chains in which the tSZ 
and kSZ terms have been fixed to their predicted, al- 
though cosmologically-rescaled, values. All four com- 
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Fig. 6. — ID likelihood curves for the tSZ power spectra in fiK^ at 150 GHz and £ = 3000 {Left panel) and ratio of measured to 
predicted tSZ power, Atgz (Right panel). In the right panel, the vertical dotted line marks the expectation value for Atsz (i0.04 due to 
sample variance). All curves assume the CSF kSZ template, and modified BB CIB frequency scaling - the tSZ results are unaffected by 
these choices. The black solid line shows the results with the Shaw tSZ model. Results for the Sehgal tSZ model are shown by the blue 
dashed line. The power constraints (left) are identical, leading to much lower A^sz's (right) as the Sehgal model predicts more tSZ power 
for a given cosmology. Finally, if we allow for possible correlations between the CIB and tSZ power spectra, the tSZ likelihood expands 
towards lower tSZ power (green line). In the correlation case, we have assumed the Shaw tSZ model - the impact is similar for the Sehgal 
model. 

hgal model. The X-ray constraint, dg ~ 0.8, which is 
independent of the tSZ modeling, rules out this line of 
argument. 



binations of the Sehgal and Shaw tSZ models and the 
NR and CSF kSZ models are considered. We show the 
change in the best-fit for each case relative to our 
baseline model (Shaw tSZ, CSF kSZ). Models that pre- 
dict more SZ power are a slightly worse fit to the data. 
Table |4] shows that the measured tSZ power is a factor 
of 0.41 ± 0.13 times that predicted by the Sehgal tem- 
plate. This template predicts a tSZ amplitude that is 
poorly matched to the data, resulting in a Ax^ — +8.4 
relative to the baseline model. Likewise, the NR kSZ 
template predicts 1 fjK^ more power than the CSF kSZ 
template and thus provides a worse fit to the measured 
powers, increasing by 2.1 (for Shaw tSZ) and 2.3 (Se- 
hgal tSZ). Adding a kSZ contribution from patchy reion- 
ization would further increase the model tension. 

The consistency between the SZ models and measured 
bandpowers can be tested further by including additional 
external cosmological constraints. Both the tSZ and kSZ 
power spectra scale most sensitively with ag. We there- 
fore importance sample each of the MCMC chains pre- 
sented in Table [5] imposing a prior of erg (JIm/O. 25)" '*^ ~ 
0.813±0.027, based on the X-ray clust er abundance mea- 
surements of [Vikhlinin et ah (2009). The subsequent 
change in i s given in the third c olumn of the table. 
The prior from Vikhlinin et al. ( 2009 ) increases the signif- 



7.2. (Tg constraints 

To obtain joint constraints on trg from both the pri- 
mary CMB anisotropy and the tSZ power spectrum, we 
must determine the probability of observing the mea- 
sured value of j4tsz for a given set of cosmological param- 
eters. In the absence of sample variance ~ and assuming 
the predicted tSZ template to be perfectly accurate - 
the constraint on ag would be given by its distribution 
at Atsz = 1 for each of the models. In practice, the ex- 
pected sample variance of th e tSZ signal is non-n egligible. 
Based on the simulations of Shaw et al. (2009), LIO as- 
sumed the sample variance at I = 3U00 to be a lognor- 



mal distribution of mean and width 



(Atsz) 



0.12. 



icancc of the discrepancy with models that predict more 
tSZ power. This result is particularly interesting since it 
eliminates one possible solution for the high-power tSZ 
models - massive neutrinos. The CMB data alone is con- 
sistent with high n eutr ino masses which slow the growth 
of structure (see [7.3). If the sum of neutrino masses 
was around 1 eV and as — 0.6, the measured tSZ power 
would in fact be higher than that predicted by the Se- 



The survey area analyzed here is eight times that of LIO 
and thus we reduce the width of this distribution by \/8- 
For each of the tSZ templates, we construct new MCMC 
chains by importance sampling to include this prior. 

The ID (Tg constraints for the two tSZ templates are 
shown in Figure [7] and Table |4j Using the Shaw template 
results in erg — 0.797 ± 0.011, while the Sehgal model 
prefers cg = 0.767 ± 0.009. The constraint without the 
tSZ amplitude information (i.e. from the primary CMB 
+ BAG + Ho alone) is erg = 0.812 ± 0.018. Adding the 
tSZ information therefore both reduces the size of the 
statistical error by 40% and lowers the preferred value, 
depending on the amplitude of the tSZ template. 

The difference in the preferred value of erg provided 
by the Shaw and Sehgal templates is significantly larger 
than the statistical errors. This implies that the system- 
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TABLE 4 

Atsz, A^SZ CONSTRAINTS 



Model 


^tsz 


In Asz 


^kSZ 


a"8 [Atsz — 1] 


Shaw, CSF, modified BB 
Sehgal, CSF, modified BB 
Shaw, CSF, modified BB + tSZ-CIB 


0.70 ±0.21 
0.41 ± 0.13 
0.60 ±0.24 


-0.37 ±0.30 
-0.89 ±0.31 
-0.54 ±0.41 


< 1.68 

< 1.57 

< 3.98 


0.797 ±0.011 
0.767 ±0.009 
0.796 ±0.012 



Note. — Atsz, Aj^sz, and erg constraints. The distributions shown in Figme[6]are clearly non-Gaussian, therefore we give 
constraints on logAtsz as well. Constraints on Aksz are 2 a upper limits. Constraints on erg include a prior on logAtsz ~ 
0.00 ± 0.04; the width of this prior reflects sample variance but not model uncertainty. The CMB±BAO±Ho constraint on as 
is 0.812 ±0.018. 



TABLE 5 
Delta for SZ models 



tSZ Model 


homog. kSZ Model 




Ax' 








w. X-ray 


Shaw 


CSF 


(0) 


0.0 


Shaw 


NR 


2.1 


2.3 


Sehgal 


CSF 


8.4 


13.0 


Sehgal 


NR 


10.7 


16.2 


Free Shaw 


Free CSF 


-1.9 


-1.9 



Note. — Change to the best-fit when the tSZ and kSZ 
contributions are fixed to the predictions of different models. 
Differences are reported relative to the Shaw tSZ / CSF kSZ 
case. The model predictions are scaled to account for the 
speci fic cos mology at each chain step as described in §5.2| 
and §5.3. 1| Note that the last row has an additional two 
free parameters describing the amplitudes of the tSZ and kSZ 
power spectra. 



atic error due to the theoretical uncertainty in the tSZ 
modeling significantly exceeds the statistical error. To 
produce a constraint on as that accounts for this, we 
must add an additional theory uncertainty to the sample 
variance when setting the prior on Atsz- 

It is difficult to determine the level of theory uncer- 
tainty on the predicted tSZ amplitude for a fixed set of 
cosmological parameters. A significant fraction of the 
total tSZ power is contributed by groups and clusters of 
mass Mson < 2 x lO^'^h-^ Mp, and redshift z > 0.65 



( Trac et al.||2011 
2U10p . lb date, sue 



Battaglia et al. 
1 objects have no 



2011b Shaw et al. 



been studied in de- 



tail with SZ and X-ray observations. From a theoretical 
perspective, models and simulations have demonstrated 
that astrophysical processes such as feedback from ac- 
tive galactic nuclei and turbulence-driven non-thermal 
pressure support can significantly alter predicted clus- 
ter thermal pr essure pr ofiles, and t hus t heir SZ signal 



2011 Parrish et al. 



(pattaglia et al. 2011a; Trac et al 
|2bll . .Shaw et al. 2010) . 

Fmlowmg yil, we choose to add a 50% theory uncer- 
tainty on In(Atsz)- This uncertainty encompasses all four 
models considered by Sll. As the sample variance is in- 
significant compared to the theory uncertainty, the new 
prior on Atsz is effectively a lognormal distribution of 
width cr(ln ^tsz) = 0.5. With this assumption, we obtain 
constraints of as = 0.807 ± 0.016 for the Shaw model. 

We estimate the contribution of the systematic the- 
ory uncertainty to the total uncertainty quoted above 
as follows. We form a new prior based on the median 
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Fig. 7. — erg constraints for the Shaw (black solid line) and 
Sehgal (blue dashed line) tSZ templates. No model uncertainty 
has been included. 

value of Atsz = 0.7 with the sample- variance based width 
of o'(ln(Atsz)) = 0.042. By sampling around the pre- 
ferred value of Atsz instead of unity, we avoid introduc- 
ing artificial tension if the modeled amplitude is incor- 
rect. Importance sampling the MCMC chain with this 
prior leads to parameter constraints that include only 
statistical noise and sample variance. The resulting as 
constraint is a (as ) = 0.010. Assuming the errors add 
in quadrature, this result implies a systematic modeling 
uncertainty of a{as) = 0.012. 

Given the current modeling uncertainties, the mea- 
sured amplitude of tSZ power in this work only 
marginally improves our knowledge of as over existing 
constraints. A detailed and systematic study of the ther- 
mal pressure profiles of low-mass (« 10^*/i~^ -^o) clus- 
ters at 2 « 0.65 is required in order to reduce the theory 
uncertainty and thus increase the precision of as con- 
straints derived from the tSZ power spectrum. 

7.3. constraints 

We can apply the methodology outlined above to other 
cosmological parameters, most notably the sum of the 
neutrino masses, ^to^. For CMB ± Hq ± BAO, neu- 
trino masses are highly degenerate with as as shown in 
Figure [8] Higher neutrino masses lead to lower as since 
massive neutrinos slow the growth of structure below the 
neutrino free-streaming length. Introducing massive neu- 
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trinos weakens the CMB + Hq + BAO constraint on tJs 
to CTg = 0.756 ± 0.044. The sum of the neutrino mass is 
constrained to be less than 0.52 eV at 95% CL. 

The tSZ power spectrum presents an independent 
probe of tig that breaks the ^ m^, degeneracy, and 
thereby improves the neutrino mass determination. The 
preference for positive neutrino masses increases with 
models that predict more tSZ power. Models that pre- 
dict more tSZ power require a lower erg to match the 
observed tSZ power. As can be seen in Figure |8) lower 
CTg values favor larger neutrino masses. 

With massive neutrinos and no modeling uncertainty 
(but including a 4% sample variance), adding the tSZ 
information reduces the uncertainty on erg by a factor 
of 2 - 3. The median erg moves up slightly to cg = 
0.776 ± 0.019 with the Shaw model and down slightly 
to CTg = 0.732 ± 0.017 with the Schgal model. In this 
optimistic scenario, we would find ^ — 0.15 ± 0.09 
or 0.29±0.10eV for the Shaw and Sehgal models, respec- 
tively. Note that the X-ray CTg measurement mentioned 
earlier is in tension with the preferred CTg with the Sehgal 
model. With a more realistic 50% modeling uncertainty 
on the Shaw model, the constraint is CTg = 0.768 ±0.031. 
With this modeling uncertainty, the Sehgal and Shaw 
model results are consistent at ~ 0.5 ct. The upper limit 
on J^^u is reduced from 0.52 eV to 0.40 eV at 95% CL 
with the addition of tSZ information. Even with a large 
theory prior, the SPT tSZ measurements improve the 
neutrino mass constraints by disfavoring regions of pa- 
rameter space with very low CTg values. 

We consider how the uncertainties on CTg and J2 
break down between systematic and stati stical uncer- 
tainties using the method presented in \ 7.2 The prior is 
constructed as before, and centered the preferred value of 
^tsz = 1.14 in the ^ rriu > eV chain. We find that the 
uncertainty on CTg can be allocated as CTsys(CTg) = 0.024 
(for the modeling uncertainty) and CTstat(o'g) = 0.020 (for 
statistical noise plus sample variance) . Similarly, the un- 
certainty on ^ TOj/ can be divided into CTsys(X!"^i') ~ 
O.lOeV and CTstat(ETO^) = O.lOeV. The 50% theory un- 
certainty is dominant for both parameters. 

8. KSZ INTERPRETATION 

We next turn to comparing our measurements of the 
kSZ power spectrum with the predictions of recent mod- 
els and simulations. We first compare our constraints 
on the amplitude of the kSZ power spectrum with our 
homogeneous kSZ templates. We then briefly mention 
the implications of the kSZ measurement for the epoch 
of reionization. Zll discuss what we can learn about the 
epoch of reionization in more detail. 

As previously presented, in the baseline model, we find 
the kSZ power to be robust to template shape between 
the CSF, NR, and patchy kSZ templates. In contrast 
to the tSZ, we see a small dependence on the frequency 
scaling of the CIB with slightly lower kSZ upper limits 
for the Single SED scaling than for the baseline modified 
BB model. However, the major modeling uncertainty's 
with the tSZ, is in the tSZ-CIB correlation (see FigureM. 
This correlation is highly degenerate with the kSZ power 
(see Figure [5]) and, when introduced as a free parameter, 
nearly triples the allowed kSZ 95% upper limit. The kSZ 
template shape affects the kSZ upper limit at the ^ 20% 
level when tSZ-CIB correlations are allowed. 



8.1. Model comparison 

In an identical manner to the tSZ, we define a dimen- 
sionless scaling parameter A^gz for the homogeneous kSZ 
power . We use the NR and CSF kSZ models of |Shaw| 
et al. (20111. Equation 13 gives the scaling of the CSF 



kSZ template $|ooq with cosmological parameters. 

In Figure |10[ we show the constraints obtained on 
^ksz- The solid black line shows our baseline model 
using the Shaw tSZ template, the CSF kSZ template 
and the modified BB CIB frequency scaling. The blue 
dashed line shows those obtained when a free tSZ-CIB 
correlation is included in our modeling. The numerical 
constraints are also given in Table |4j Recall that a value 
of ^ksz — 1 implies that the measured kSZ power is con- 
sistent with the kSZ template predictions, while Atsz less 
(greater) than 1 indicates that the model overestimates 
(underestimates) the amplitude. Note that sample vari- 
ance will cause deviations away from unity on the order 
of 0.01. 

In this work, we have significantly tightened the up- 
per limits on the kinetic SZ power. These limits are 
not in significant tension with o ur fiducial homoge neous 
kSZ template (the CSF model of [Shaw et aL][201l| . The 
2ct upper limit on Aksz exceeds 1 tor all of the MCMC 
chains in Table |4j However, for our baseline model, the 
upper limit is only modestly in excess of 1. The Single 
SED CIB model allows for slightly less kSZ power than 
the modified BB frequency scaling. Therefore, while the 
data are consistent with the CSF kSZ prediction, there 
is little room for an additional contribution from patchy 
reionization. Zll note that, if the duration of reioniza- 
tion spans Az = 2 we expect an additional 1 of 
kSZ power. For Az = 4 we expect 2/iK^. Clearly when 
combined with the post-reionization signal, the predicted 
kSZ will significantly exceed that measured by SPT in 
the baseline model. 

incorporating a tSZ- 



However, as discussed in j 6.2.1 



CIB spatial correlation into our model allows for consid- 
erably more kSZ power than our fiducial model. In this 
case, the upper limit on A^sz is 4.13, which allows for 
substantially more power than predicted by the homoge- 
neous kSZ models alone, and thus a significant contribu- 
tion from patchy reionization. 

8.2. Patchy reionization 

We now switch to considering limits on the kSZ sig- 
nal contributed by patchy reionization. For this we use 
a patchy kSZ template from Zll based on a model in 
which reionization started at z = 11 and ended z = 8. 
As would be expected, the kSZ signals from the reion- 
ization and post-reionization era are nearly perfectly de- 
generate. Therefore, our constraints on the patchy kSZ 
power depend on the modeling of the post-reionization 
homogeneous kSZ power. We choose to bracket the latter 
with three cases: no homogeneous kSZ, the intermediate 
CSF homogeneous kSZ model, and the highest NR ho- 
mogeneous kSZ model. 

The constraints on patchy kSZ power are shown in 
Figure |9] Note that we see less difference in the al- 
lowed patchy kSZ power than would be inferred from 
the power differences between the assumed homogeneous 
kSZ model. The smaller than expected differences are 
because the likelihood function for the total kSZ power. 
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Fig. 8. — Left panel: 2D likelihoods for erg and J^m^. The filled colored contours show the 1, 2, and Scr constraints for the CMB 
+ HO + BAO data. The black solid lines show the constraints with the addition of the tSZ information (assuming a 50% theory prior). 
Right panel: ID likelihoods for the sum of neutrino masses for four cases. The red dot-dash line shows the constraints without including 
tSZ data. The black solid line shows the constraints when we include the measured tSZ power with the Shaw tSZ model and a 50% 
theory prior. This tightens the upper limit on the sum of the neutrino masses from 0.52 eV to 0.40 eV at 95% CL. We show the constraints 
on neutrino mass that would be derived with perfect tSZ modeling in the final two lines: with the Shaw tSZ model (purple dashed line) 
and Sehgal tSZ model (green dotted line). The Sehgal model predicts more tSZ power for a given erg. As a result it pushes the allowed 
parameter volume to lower erg and higher ^ m^, leading to a nominal detection of the total neutrino mass at 0.29 it O.lOeV. 



before any amount of homogeneous kSZ power is sub- 
tracted, already peaks at zero. We find 95% CL upper 
limits on the kSZ power due to patchy reionization of 2.6, 
2.0, and 1.9 /iK^ for the zero, CSF, and NR homogeneous 
kSZ models respectively. As with all the kSZ constraints, 
we find the patchy kSZ constraints are most sensitive to 
tSZ-CIB correlations. With a free tSZ-CIB correlation 
and the same set of three homogeneous kSZ models, the 
95% CL upper limits on patchy kSZ power are degraded 
to 5.7, 4.9, and 4.7 /iK^. Zll interpret these upper limits 
on patchy kSZ power to constrain the duration of the 
epoch of reionization and the ionization history of the 
Universe. 

9. COSMIC INFRARED BACKGROUND CONSTRAINTS 

The combination of Planck/KFl and SPT data 
strongly constrains the CIB power spectrum at millime- 
ter wavelengths, especially at 220 GHz. The CIB con- 
straints for a variety of model assumptions are sum- 
marized in Table [6] and Figure 11 

cases, sources with flux >6.4mJy at 150 GHz have been 
masked. We find the 220 GHz Poisson amplitude at 
i = 3000 to be DPqqo = 68.0 ± 3.3 ^K^. The cluster- 
ing amplitude is D^qqq = 56.3 ± 4.2 /iK^ - about 45% 
of the total CIB power at this angular scale. The cross- 
over point where the Poisson and clustering terms have 
equal power is at £ = 2500. At 150 GHz, the Poisson 
and clustering powers are -Dfooo — '^■^'^ ^ 0.38 /l(K^ and 



ZJfooo = 6.25 ±0.52 /xK^ respectively. The estimated GIB 
at 95 GHz is very small, but this is an extrapolation from 
higher frequencies and sensitive to the frequency model- 
ing. Both Poisson and clustered CIB components have 
more power at 150 GHz than the SZ spectra across the 
multipoles {£ > 2000) measured in this work. We will 
leverage the wide range of angular scales probed by SPT 
and Plan ck/H Fl to constrain the shape of the GIB clus- 
tering in §9.1[ 

Both models considered for the frequency scaling of 
the CIB lead to similar results. The data require a steep 
drop-off in CIB power from 220 to 150 GHz; the spectral 



index is a 



150/220 



= 3.51 ±0.04 and 



ai50/220 



= 3.56±0.07 



for the Single SED and modified BB models respectively. 



This is shallower by 1 cr than the ai 



50/220 



3.68 ±0.07 



measured fo r the clustered spectral index by Addison] 
et al I (|2011[ ) using data from BLAST, Planck/RFl, and 
AC'l'. We note that when we allow the Poisson and clus- 
tered spectral indices to differ later in this section, the 
inferred clust ered spectral index b ecomes more consis- 



Recall that, in all ^gj^^- ^j^i^ thelAddison et al.lpUni) result. The observed 



spectral index tor the lowest two published Planck/HFl 
channels (217 and 353 GHz) is shallower: 0217/353 — 
3.22 ±0.08. 

These indices are steeper than many current models 
predict. For instanc e, the best-fit parameters found by 
when fitting COBE F IR AS data 
2.6. The IPlanck Collaboration 



Fixsen et 
predict 



al. (1998) 



0^150/220 
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et al. (20111 previously pointed out that these param- 
eters also overpredict the CIB anisotropy fro m 217 - 
857 GHz. The alternate parameter fitting by 'Gispert' 



et al. ( 2000) that works well across the Pl anck/HFl ire- 
quencies ( Planck Collaboration et al.|2011[ ) still predicts 
a shallower q;i5o/ 22o = '^■04 between 150 and 220 GHz. 



Fig. 9. — ID likelihood curves for the kSZ power (in /^K^). Left panel: Constraints on the kSZ power with different CIB models and 
kSZ templates. The modified BB CIB frequency scaling and Shaw tSZ template are assumed in all cases. Solid lines denote constraints 
without tSZ-CIB correlations; dashed lines mark constraints with tSZ-CIB correlations allowed. Black lines show the results with the 
CSF kSZ template (the NR kSZ template is similar). Red lines mark the constraints when we alternatively use a template based on the 
kSZ contribution from the epoch of reionization. Right panel: The constraints on the patchy kSZ power when we include a model 
prediction for the homogeneous kSZ power. Solid lines denote constraints without tSZ-CIB correlations; dashed lines mark constraints 
with tSZ-CIB correlations allowed. Black lines show the results when the CSF homogeneous kSZ model is assumed. This model includes 
cooling and star formation which suppresses the kSZ power. Red lines mark the constraints when the higher amplitude NR homogeneous 
kSZ model is assumed. 

data - a2i7/353 = 3.22 ±0.08. The changing slope means 
that the modified BB model prefers lower temperatures 
(16 ± 6K, with a non-Gaussian distribution) and high 
values of /? (/3 = 1.80 ± 0.14) when we include the 353 
GHz data. The modified BB parameters, with and with- 
out 353 GHz data, are consistent at < la. 

With the Single SED frequency scaling, the median 
Zc = 0.53, but the best-fit Zc ~ 0. Recall that in the 
Single SED model, the DSFG density as a function of 
redshift is described by a Gaussian centered at z^. The 
distribution of z^ tightens around zero when the 353 GHz 
data are added. This suggests that either the adopted 
SED template is not representative or that using a sin- 
gle SED for all galaxies is an over-simplification. Note 
that in the Single SED model, the maximum allowed 
spectral index is a = 3.57 at Zc = 0. For the specified 
model parameters, higher spectral indices would require 
Zc < which is not allowed by the z^ prior. The data 
are exploring the area at this limit, which slightly low- 
ers the median spectral index for the Single SED model. 
These redshifts are driven low to accommodate the steep 
observed spectral index. These constraints are shown 



Bethermin et al. (2011) modeled DSFG source counts 
using data from IKAH, Spitzer, Herschel, and AzTEC 
- these modeled source counts were used to predict 
the Poisson power in the P lanck/HFl CIB bandpowers 



(Planck Collaboration et al.,2011). The spectral index 



predicted by the Bethermm model is 01150/220 = 2-8. The 
maximum spectral index among the ten provided realiza- 



tions of the Bethermin model is a 



150/220 



3.1. These 



models are a poor fit to the observed CIB power at mil- 
limeter wavelengths. For instance, the Bethermin model 
predicts ~25% too much Poisson power at 220 GHz and 
twice the observed power at 150 GHz. 

There is a clear degeneracy between the two param- 
eters of the modified BB model with only < 220 GHz 
data. Recall that in this model, the CIB spectrum is 
described by i^/^B^(T). Without 353 GHz data, temper- 
atures above ~ 15 K are constrained only by the prior 
since the observed frequencies are in the RJ tail of the 
spectrum. The preferred temperature is > 12 K, and 
/3 = 1.8 ± 0.1. These results are plotted in Figure [12} 
Adding 353 GHz data partl y bre aks the T-f3 degener- 
acy. We showed earlier in { 5.4.3| that the spectral in- 
dex was shallower in the higher frequency Planck/HFl 



in Figure 12] The tighter spectral index range with the 
Single SED model also (artificially) tightens the kSZ con- 
straints by a small amount. 

We explore scaling the clustered and Poisson CIB in- 
dependently with frequency, and find a small prefer- 
ence for different spectral indices when considering only 
< 220 GHz data. For instance, with the modified BB 
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^kSZ 

Fig. 10. — Poste rior distribution of the kSZ scaling parameter 
^kSZ (see Equation [20] and Table [ijl for different CIB models. The 
vertical dotted line marks the model expectation value for Aj^gz 
(±0.01 due to sample variance). The CSF kSZ model we consider 
here is the lower of the two kSZ models and does not include the 
kSZ contribution from patchy reionization. Therefore, Ajjgz > 1 
would be reasonable. The solid black line shows the baseline 
model using the Shaw tSZ template, the CSF kSZ template, and 
the modified BB CIB frequency scaling. The blue dashed line 
shows the results obtained when the tSZ-CIB correlation is allowed 
in our modeling. 

model, the 150 - 220 GHz spectral index is ai5o/220 = 
3.45 ± 0.11 and 3.72 ± 0.12 for the Poisson and clus- 
tered term respectively - a difference of 1.7 a. The spec- 
tral index constraint is modestly weakened by the addi- 
tional free parameter, but both parameters remain well- 
constrained. The SZ constraints are not affected. How- 
ever, we note that adding the 353 GHz Planck/HFl data 
removes this small preference for different spectral in- 
dices. In this case, the Poisson and clustered spectral 
indices are 0159/220 = 3.46 ±0.09 and 3.50 ±0.08 respec- 
tively. 

9.1. Angular dependence of the clustered CIB 
component 

As noted earlier, we find strong evidence for clustered 
DSFGs in the data. We have modeled the clustering by 
a power law with cx i'^'^; here we examine if the data 
prefer an alternate shape. The SPT data alone can be 
fit with either a linear-theory clustering template or the 
power-law baseline template - only the extremely wide 
angular range probed by the combination of SPT and 
HFI allows shape discrimination in the clustered tem- 
plate. We parametrize the CIB clustering shape in two 
ways: fitting for the exponent of the power law, or fitting 
for the amplitude of a linear-theory term in addition to 
the power law term (with fixed exponent). Adding lin- 
ear theory power is similar in effect to lowering the power 
law exponent. 

We introduce the power law exponent, 7, as a new 
parameter to the MCMC parameter chains. A uniform 
prior is imposed from 0.1 to 1.4. The upper edge is chosen 
to prevent overlap with the l'^ scaling of the Poisson term. 
The constraints on 7 are shown in Figure 13 The data 



prefer 7 = 0.81 ± 0.08, consistent with the nominal value 
of 0.8. Adding the Planck/RFl 353 GHz CIB bandpowers 
improves the shape measurement to 7 = 0.82 ±0.05. The 
likelihood of the best-fit point is not improved by the 
introduction of this new parameter for either set of data. 
Both results are consistent with the recent constraint of 



7 = 0.75 ± 0.06 b y [Addison et al.| (|2011[ ) using data from 
BLAST, Planck/iiFl, and AC'l'. For simplicity and ease 
of comparison with earlier results, we continue to use 
7 = 0.8 in all other parameter chains. 

We also test introducing the linear theory model with 
a free amplitude in addition to the power law term with 
7 = 0.8. The constraints on the amplitudes of the linear 
theory and power law terms at 220 GHz and i = 3000 
are shown in Figure 14 We do not find a significant 



preference for a non-zero amplitude of the linear the- 
ory term; the derived 1 a confidence intervals are [-0.02, 
0.08 /iK2], [-0.2, 0.7^iK% and [-1.4, 6.2^K2] at 95, 150, 
and 220 GHz respectively. We have allowed negative val- 
ues to allow the maximum shape flexibility; the total 
power in the two clustering terms is always positive. Note 
that these ranges should be compared to 0.67 ± 0.13, 
5.7 ± 1.0, and 51 ± 9 /iK^ for the corresponding ampli- 
tude of the power-law term. The quality of the flt is not 
improved by adding this free parameter. The preferred 
values imply that most clustered CIB power is in the 
power law shape instead of the linear theory shape at 
£ > 500. Adding the Planck/BFl 353 GHz bandpowers 
does not significantly change this picture. 

9.2. CIB constraints with tSZ-CIB correlation 

The galaxy clusters that contribute to the tSZ power 
spectrum have a higher than average density of galax- 
ies, each with some amount of star formation. This will 
lead to an anti-correlation between the tSZ and CIB 
power spectra. W c intro duce this correlation following 
the treatment in j 5.4.4 Note that in this treatment, 
the correlation is defined with respect to all CIB power, 
although the Poisson and clustered components may cor- 
relate differently. As discussed in Zll, some models pre- 
dict a correlation effectively in dependent of ^ at current 
instrumental sensitivities (e.g. Sehgal et al.||2010 ), while 
others predict that the magnitude of the correlation will 
rise across the relevant multipoles. We use a single cor- 
relation coefficient on all angular scales. The data prefer 
a correlation coefficient of ^ = — 0.18±0.12 (i.e., an anti- 
correlation between tSZ and CIB signals, consistent with 
a fraction of the CIB emission being spatially coincident 
with the tSZ decrement from clusters). This constraint 
depends on the assumption that the correlation is inde- 
pendent of £. Wc tested a correlation model described 
by Zll where the absolute magnitude of the correlation 
doubles from i = 2000 — 10000, and found similar uncer- 
tainties but a mean value 1 a closer to zero a,t i — 3000. 
In the near future, we expect observations with Herschel 
and Planck in combination with SPT to place better ob- 
servational constraints on this correlation. 

Introducing a tSZ-CIB correlation does very little to 
the inferred CIB powers; however, the uncertainties in- 
crease by 10-50%. With tSZ-CIB correlations, the con- 
straint on the Poisson power at 150 GHz is D^qqq = 
8.04±0.48 /iK^, and the constraint on the clustered power 



at 150 GHz is D 



3000 



6.71±0.74/xK^ The best-fit spec- 
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TABLE 6 
CIB Constraints 



Model 


95 GHz 


Poisson (/.tK^^ 
150 GHz 


220 GHz 


150 — 220 


^150-220 tSZ-CIB correlation 


Mod. BB 


0.87 ±0.09 


7.54 ±0.38 


67.96 ±3.29 


3.56 ±0.07 


- 




1.02 ±0.15 


7.97 ±0.49 


65.99 ±3.65 


3.45 ±0.11 


3.72 ±0.12 


Linear th. clustering 


0.90 ±0.10 


7.73 ± 0.48 


69.10 ±3.88 


3.55 ±0.08 


- 


Power law 


0.86 ±0.10 


7.50 ±0.52 


67.79 ±4.23 


3.57 ±0.08 


- 


tSZ-CIB corr. 


1.04 ±0.15 


8.05 ±0.48 


66.67 ±3.68 


3.45 ±0.11 


-0.18 ±0.11 


Single Sed 


0.87 ±0.05 


7.74 ±0.34 


66.51 ±2.96 


3.51 ±0.04 


- 




Non-linear Clustering (/xK"^) 


Lin. 


Theory Clustering (/^K"^) 




95 GHz 


150 GHz 


220 GHz 


95 GHz 


150 GHz 220 GHz 


Mod. BB 


0.72 ±0.09 


6.25 ±0.52 


56.26 ±4.22 






O-c 7^ Op 


0.59 ±0.11 


5.77 ±0.59 


57.26 ± 4.34 






Linear th. clustering 


0.67 ±0.13 


5.73 ± 1.00 


51.13 ±9.12 


0.03 ±0.05 


0.26 ±0.43 2.39 ±3.78 


Power law 


0.72 ±0.11 


6.28 ±0.71 


56.70 ±6.39 






tSZ-CIB corr. 


0.87 ±0.17 


6.73 ±0.73 


55.85 ±4.47 






Single Sed 


0.73 ±0.06 


6.46 ± 0.48 


55.48 ±4.10 







Note. — Constraints on the CIB power and spectral index are presented for six CIB models. Five of these models use the 
modified BB frequency dependence for the CIB. The first row is for the baseline model. The next four rows consider 1-parameter 
extensions to the baseline CIB treatment. In the "ac 7^ Qp" row, the Poisson and clustered terms have independent frequency 
scalings. In the "Linear th. clustering" row, we allow for a (positive or negative) linear theory clustering term in addition 
to the power law clustering term. In the "Power law" row, we introduce the exponent of the power law clustering term as 
a free parameter instead of fixing it to 0.8. In the "tSZ-CIB corr." row, we allow for correlations between the tSZ and CIB 
components. The last row considers constraints under the alternate Single SED frequency modeling for the CIB. For each model, 
we present the constraints on the amplitudes in /iK^ at £ — 3000 of the Poisson and power-law clustering CIB templates. For 
the linear-theory clustering row, we also present the preferred amplitudes of the linear-theory term - none of which represent 
detections. Finally, we present the effective spectral index from 150 to 220 GHz for the Poisson (and if different) the clustering 
terms. 



tral index is lower by 1 a, and the uncertainty increases 
by 50% to give a = 3.45 ± 0.11. However as discussed 
earlier, a free correlation significantly degrades the kSZ 
constraint. 

10. CONCLUSIONS 

We have presented the CMB temperature anisotropy 
power spectrum from the complete 800 deg^ of sky ob- 
served by SPT in the 2008 and 2009 austral winters. The 
six cross-spectra from observations at three frequencies 
(95, 150, and 220 GHz) are shown in Table|l| The band- 
powers measure arcminute-scale {£ > 200U) millimeter 
wavelength anisotropy. On these small scales, DSFGs, 
radio galaxies, and the kinetic and thermal SZ effects 
contribute significant power in addition to the primary 
CMB anisotropy. 

Using MCMC methods, we perform multi-frequency 
fits to a combined data set including WMAP7, BAO, 
Hq, previously published low— £ SPT bandpowers, 
Planck/RFl CIB bandpowers, and the e > 2000 SPT 
bandpowers presented in this work. We find that the 
minimal model needed to explain the current data adds 
five free parameters beyond the six ACDM parameters: 
four describing DSFGs and one describing the amplitude 
of the tSZ power spectrum. We also include a free pa- 
rameter for the amplitude of the kSZ spectrum, and a 
prior on a Poisson population of radio galaxies. This 
six parameter extension improves the likelihood of the 
best-fit model by Ax^ = —2081. We explore a number 
of extensions to this baseline model, including different 
templates for the kSZ, tSZ, or CIB terms, different mod- 
eling of the CIB frequency spectrum, and correlations 



between the tSZ and CIB. We do not find these exten- 
sions significantly improve the quality of fits to the data, 
however, they may slightly alter the conclusions drawn. 

We fit for the amplitude of the tSZ and kSZ power 
spectra, using a single amplitude parameter defined in 
units of ^K^ at i 3000 and 150 GHz. The constraint 
on the 150 GHz tSZ power is ZJ^ggp = 3.65 ± 0.69 /iK^ 
in the baseline model. We find the tSZ amplitude is 
independent of the template shape for the two models 
considered: the Shaw and Sehgal tSZ models. We find 
the kSZ results are also robust with respect to template 
shape. The 95% confidence upper limit on the kSZ power 
is i?|sz^<2.8/iK2. 

We test whether extensions to the primary CMB mod- 
eling, including running of the scalar index, changing 
number of neutrino species, or massive neutrinos, af- 
fect the SZ constraints and find they do not. We also 
investigate whether alternative CIB modeling assump- 
tions lead to changes in the SZ constraints. The SZ 
constraints are fairly insensitive to the CIB modeling. 
An exception is their sensitivity to the modeling of the 
tSZ-CIB correlation. When allowing for the possibil- 
ity of this correlation, the allowed parameter space ex- 
pands along the degeneracy line of D3000 + ^■'^^sooo- 
The intersection of this line moves by less than 1 a: 
Dfooo + 0-5-D|^o^o = 4.15 ± 0.56 /.tK^ without correla- 
tions and 4.60 ±0.63 fiK^ with free correlations. We note 
that this degeneracy line closely resembles that found by 
LIO, who used a linear combination of the 150 GHz and 
220 GHz maps to subtract the CIB. When correlations 
are allowed, the 95% confidence upper limits on the to- 
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Fig. 11. — ID likelihood functions for the powers in fiK^ of the two CIB components at £ = 3000 in the three SPT bands. The top row 
shows the constraints on Poisson power; the bottom row shows constraints on the power-law clustering power. Sources detected with flux 
> 6.4 mjy at 150 GHz are not included. From left to right, the columns display the constraints at 95, 150, and 220 GHz. The likelihood 
functions with the fiducial, modified BB CIB frequency scaling are shown by the black lines. The green lines display the likelihood 
functions if the clustered and Poisson DSFG components are allowed to scale independently with frequency. If we use the Single SED CIB 
frequency scaling instead, we find the blue dashed lines. The CIB constraints for the Single SED and modified BB scalings are very 
similar. The CIB powers, particularly the extrapolation to 95 GHz, change somewhat when the Poisson and clustered component frequency 
scalings vary independently. The preferred spectral index is slightly shallower for the Poisson than clustered component. The Poisson and 
power-law clustering terms are comparable on angular scales near £ = 2500. 



tal kSZ power are D^^^^ < 6.7/^X2 and D^^^^ < 5.7 fiK^ 
for the homogeneous and patchy reionization templates 
respectively. The real kSZ signal is expected to be an 
admixture of the two templates. 

We use the measured tSZ power to constrain erg and 
the sum of the neutrino masses. For massless neutri- 
nos and the optimistic assumption of assuming perfect 
tSZ modeling, the tSZ power constraint tightens the 
CMB-fBAO-l-Ho as constraint by 40%. However, the 
preferred erg value is 3 cr discrepant between the two 
(Shaw and Sehgal) tSZ models considered. With a con- 
servative 50% modeling uncertainty, adding the tSZ in- 
formation does not significantly improve the constraints 
on (Tg . The impact of the tSZ power measurement is more 
significant when neutrinos are allowed to have mass. 
With massive neutrinos, even with a 50% modeling un- 
certainty, the tSZ power spectrum reduces the erg uncer- 
tainty by 30% and tightens the upper limit on the total 
neutrino mass from 0.52 eV to 0.40 eV. Precise measure- 
ments of the pressure profiles of high redshift clusters are 
essential to reduce the theory uncertainty and realize the 
statistical power of the tSZ measurement. 

We also use the SPT data to constrain the astro- 
physical foregrounds (after masking the brightest point 
sources) . Two parameters - which can be fixed with pri- 
ors from other observations - fit the radio source com- 
ponent. We find a four parameter model is adequate to 



describe the Poisson and clustered DSFGs. The spec- 
trum of the DSFGs falls sharply from 220 to 150 GHz, 
with an effective spectral index of a = 3.56 ± 0.07. 
This is steeper by ~2 a than extrapolations from mod- 
els based on higher frequency observations of the CIB. 
However, it is in agreement with previous spectral in- 
dex rneasure ment s at these frequencies by Sll, Dunkley] 



et al. (2011 1, and Addison et al. (2011 1. We find Poisson 



DSFG powers at ^ = 3000 of D^^= 7.54 ± O.SSfiK"^ 



68.0 ± 3.3, /iK^ at 150 and 220 GHz respec- 
tively. We find clustered DSFG powers at ^ = 3000 of 



and Dfooo 



3000 



6.25 ± 0.52 ^K^ and D] 



looo = 56.3 ± 4.2 /zK^ at 
150 and 220 GHz respectively. We find some preference 
for different spectral indices for the clustering and Pois- 
son DSFG components when we exclude the Planck/HFl 
353 GHz data; this preference vanishes when the 353 GHz 
data is included. 

With Planck/RFl and SPT data, the CIB is con- 
strained on angular scales from £ — 80 — 9400 at 220 GHz. 
This large lever arm allows both the separation of the 
Poisson and clustered DSFG components, and a tight 
constraint on the shape of the clustering component. We 
find that the data prefer the clustered power to be pre- 
dominately in the form of a power law P rather than the 
linear theory prediction. This is strong evidence for non- 
linearities in the DSFG clustering. We further constrain 
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Fig. 12. — CIB frequency modeling. Left 'panel: Constraint on the mean source redshift with the simple Single SED model parametrized 
with T = 35 K, /3 = 2, and = 2. The redshift constraints with and without the HFI 353 GHz data included are plotted in blue dashed 
and black lines respectively. Middle panel: 2D likelihood surface for T and /3 for the modified black-body CIB spectrum with 95-220 GHz 
data (shaded contours) and 95-353 GHz data (line contours). Right panel: The effective CIB spectral index for 150 - 220 GHz. For 
the Single SED model frequency scaling, the preferred spectral indices are plotted in black solid and blue dashed lines for 95-220 GHz 
and 95-353 GHz respectively. The sharp cutoff near a = 3.55 is the Zc > prior in the Single SED model. The spectral index constraints 
with the modified black-body model with and without 353 GHz are shown by the red dot-dashed and green solid lines respectively. 
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Fig. 13. — ID likelihood function for the power-law exponent, 
7, of the CIB clustering term {Di oc £'^) . The fiducial model value 
is 0.8. The curves for constraints with and without including the 
HFI 353 GHz CIB bandpowers are plotted with dashed blue and 
black lines respectively. 
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Fig. 14. — 2D likelihood function for the amplitudes in /iK^ 
at 220 GHz of the power-law and linear-theory CIB clustering 
terms. The filled yellow to red contours show the 1, 2, and 
3 CT constraints with SPT-l-PZancfc/HFI data at 95 - 220 GHz. The 
black contours show the same constraints when including the 
Planck/HFl 353 GHz CIB bandpowers. In both cases, the data 
significantly prefer the power law clustering template, and are con- 
sistent at 1 (T with zero additional linear theory power. 

the exponent of the power law to be 7 = 0.82 zfc 0.05, 



consistent with other recent constraints (Addison et al 



We test the impact of introducing a tSZ-CIB correla- 
tion coefficient, and find the data prefer a negative cor- 
relation at 98% confidence. Note that this significance 
may vary for more sophisticated CIB models. The tSZ- 
CIB anti-correlation modestly increases the inferred CIB 
powers, while counter-intuitively reducing the tSZ power 
and increasing the allowed kSZ power. The uncertainties 
on both SZ terms increase significantly when correlations 
are allowed. External data or modeling to constrain the 
tSZ-CIB correlation would be extremely valuable to fu- 
ture kSZ constraints. 

A companion paper, Zll, uses the SPT bandpowers 
and kSZ upper limits presented here to investigate the 
epoch of rcionization. Zll constrains the duration of 
reionization from the SPT measurement of kSZ power. 
The midpoint of reionization is determined from the 
WMAP constraint on the optical depth. Together, the 
measurements yield the first constraints on the begin- 
ning, duration, and end of reionization from the CMB. 

The SPT survey of 2500 deg^ finishes in November, 
2011. The survey area, comprising 6% of the total sky, 
has been mapped to depths of approximately 40, 18, and 
70 /iK-arcmin at 95, 150, and 220 GHz respectively. We 
expect bandpower uncertainties for the full dataset to 
be reduced by a factor of 1.7 compared to this work. 
We expect from simulations that the constraints on the 
tSZ, kSZ, and tSZ-CIB correlation will improve com- 
mensurately. In addition, we have selected 100 deg^ of 
the SPT survey for deeper integration and Herschel fol- 
lowup. Herschel/ SPlRFi will observe this region in the 
winter of 2011. The combined SPT / Herschel dataset will 
be an invaluable resource in understanding the CIB and 
its correlation with the SZ power spectra. We will be able 
to use the Herschel maps to clean the CIB from the SPT 
maps of this 100 deg^ field, further improving measure- 
ments of the tSZ and kSZ power spectra. These results 
will advance our understanding of structure formation 
and the reionization of the universe. 
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